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ON THE STABILITY OF INVISCID PARALLEL 
FLOW IN HYDROMAGNETICS 


by R. K. JAIN 


Department of Mathematics, University of Delhi, Delhi, India 


Summary 


In the present paper we discuss the stability of parallel flow of conducting 
and inviscid fluid between two fixed and concentric cylinders in the presence 
of an axial uniform magnetic field. It will be shown that every flow W(r) 
is stable for infinitesimal perturbations. 


§ 1. Introduction. Recently a number of papers have appeared 
on the stability of rotational Couette flow in hydromagnetics by 
Chandrasekhar 1-3), Michael®), Velikhove7’). These authors 
have considered both the presence and absence of dissipative forces. 
But the hydromagnetic stability of parallel or nearly parallel flows 
has not, so far, been given much attention. However, Stuart 8) 
and Lock #4) have obtained the curves of neutral stability for 
symmetric parallel flows with different magnetic field configurations. 
The purpose of the present paper is to discuss similar problems in 
the absence of viscosity. We shall establish that every parallel 
flow W(r) is stable, allowing only damped perturbations. 


§ 2. Formulation of the problem. We consider the parallel flow 
of an inviscid and conducting fluid between two fixed and concentric 
cylinders of radii a and b (b> a). Let (7, 6, z) denote the cylindrical 
coordinates of a point where the common axis of the cylinders has 
been taken as the z-axis. Throughout the present paper cylindrical 
coordinates will be used. The equations of hydromagnetics are 


I 
4p 


i 
DP dee (v-V)v= — — grad p + 
ot p 


curl H x H, (1) 
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oH 1 
iene Hye eR: a 
- curl(v x H) + a (2) 
div H — 0, (3) 

divy = 0. (4) 


It is clearly seen that these equations admit the stationary solution 
>= constant 17 (say);-V —10,.0,(7)|, H — (0, Ona pao, 


where H is constant. Let the perturbed state of the system (5) be 
characterised by 


b= IT 4 &, v =; Up, 0; 4 W), OS ie lg hs 


where @, Uy, Vg, Vz, Ar, hg and hz are infinitesimal quantities and are 
chosen in the form 


Vp = u(r) eee), hy = (7) etlett+kz) 
vg = v(r) enna): ho aes (7) SS 
Vz = w(r) en kay hz = x(r) euclagcey (7) 


& = &(r) etlet+ he), 


Here & is a real positive number and c is a complex quantity. The 
first order perturbation equations are given by 


Hes 
1Kpu = — Do + ae (zkp — Dy), (8) 
Ak 
Kp? =a 
oe An _ (?) 
ipKw + upW' = —ike, (10) 
é ; 1 
iKp = ikuH + 7 — (DD* — }2) 9, (11) 
. 1 ; 
Ky = tkvH + a (DD* — k?) y, (12) 
fy te ee ne ee 
: SEED SE aaah Teg" aaah 
D*u + thw = 0, (14) 


D*p + iky = 0, (15) 
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where 
d 
i ee ee 
dy dy v 
- (16) 
DD* = D*D : pp — 2 iad S 
2 dr2 y dr 
and 
kK=kW-e. (17) 


Equations (8)—(15) are not all independent. For instance letting (11) 
operate on D* and making use of (15), we shall get (13). Elimination 
of v from (9) and (12) gives 
io H2k2K* 
(DD* — k?)y = — (4x0K _ aa) y 
p [| 
Also elimination of @, w and wu from (8), (10), (11), (14) and (15) 
gives 


(18) 


ipkK pk WwW’ 
= (DD* — k2)2q — EE — 4apK2+ (we ei (DD* —k2)p + 
oO oO v 
+ 82pKkW'D*p =0. (19) 
§3. Boundary conditions. Since the boundaries are rigid, the 
normal component of velocity vanishes: 
t= Oinat— = a;.0! (20) 


We consider the boundaries to be perfectly conducting so that the 
Fermi boundary conditions 1) are valid: 


ere die 74,0, (21) 
17 iO Pattaya, 0: (22) 

From (11), (20) and (21) we get 
(DD* — k?)p=0 at r=a,b. (23) 


Equation (18) is of second order, and hence two boundary conditions 
viz. (22), are sufficient for its solution. Equation (19) is of fourth 
order and for its solution four boundary conditions, viz. (21) and 
(23), are sufficient. 
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§ 4. Proof of the theorem. The problem has been reduced to 
finding the eigenvalues of c. For this purpose we consider (18) and 
(22) | 

Equation (18), on multiplication by ry* (y* being the complex 
conjugate of y) and integration with respect to 7, gives 


b ; b H2k2K* 
i r(|D*p|2-LR? |y|2)dr+ =| H( 4xpK — ae ) wtar=o. (24) 


Adding the complex conjugate of (24) to if, we obtain 


b b H2p2 
2[ n\D*ylt + #2 [pl8) dy = [ (420 + ar) ly? dr, (25) 


where c = c, + 2¢cg. Equation (25) gives 


“y(\Dty2 +B |yl2) dr 
ee (26) 


OO eros H2b2 
IF ‘ (4 7 Ke Lark 


The two integrals on the right hand side of (26) are positive, 
establishing thereby that cg > 0. Thus every parallel flow W(r) 
allows only damped distrubances and consequently is stable. 
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A LIQUID SCINTILLATION COINCIDENCE 
COUNTER FOR RADIOCARBON *) 


by K. J. NYGAARD 


Institutt for Teknik Fysikk, Norges Tekniske Hogskole, Trondheim, Norway 


Summary 


Measurements have been made of the radiocarbon counting performance 
of a liquid scintillation coincidence counter used in conjunction with an 
anticoincidence shield counter. When operated at room temperature, the 
system has a radiocarbon counting efficiency of 59% at a background of 
16 counts/min. The activity of ethanol made from contemporary wood has 
been determined to 13.1 + 0.6 disintegrations per minute per gram carbon. 


§ 1. Introduction. The internal sample liquid scintillation counter 
admits the sample in the liquid (or solid )state. Thus it is usually 
possible to introduce larger carbon amounts into a liquid scintil- 
lation counter than into a gas counter. This is especially important 
in the radiocarbon dating of old samples, where a counter of high 
radiocarbon concentration efficiency and low background is needed. 
In conjunction with the Trondheim Radiocarbon Dating Laboratory 
an attempt was accordingly made to utilize the simple photon 
coincidence detector system of Sigmond and Schjetne!)§) for 
a room temperature operated radiocarbon liquid scintillation 
counter. 

For review of liquid scintillation counters for radiocarbon dating 
is referred to Olsson 8). 


§ 2. Counting system. The counting system is shown schemati- 
cally in fig. 1. The 100 ml scintillator cell consists of a stainless 
steel cylinder highly polished on the inside, connected by means of 
teflon sealing rings to plane quartz end windows. The 13-stage 
venetian blind photomultiplier tubes (E.M.I. 6255) are coupled to 
the windows by means of thin layers of silicone oil. 


*) Condensed from the author’s diploma thesis, NTH 1959. 
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The photon coincidence detector system is almost identical with 
that described by Sigmond and Schjetne in their paper on tritium 
liquid scintillation counting, the most notable difference being that 
the multipliers used in the present work were somewhat more noisy. 
The main features of this system are that the photomultiplirs are 
coupled directly to the coincidence circuit (the cathode followers 
used are not really necessary), and that the pulse height selection 
is performed on the coincidence circuit output pulses. This is made 
possible by the quasilinear coincident pulse transfer properties of 
the DeBenedetti-Richings coincidence circuit. For further 
details the reader is referred to Sigmond and Schjetne. 


MULTIPLIER SHIELD 
E.H.T. E.H.T. 
SUPPLY SUPPLY 


[ ANTI-COINCIDENCE SHIELD COUNTER 
CATHODE |PHOTO- SAMPLE||PHOTO- Cat 5 
FOLLOWER | MULTIPLIER ||CELL  ||MULTIPLIER | FOLLOWER 


ANTI-COINCIDENCE SHIELD COUNTER 


COINCIDENCE 
CIRCUIT 


LINEAR 
AMPLIFIER 


LINEAR 
AMPLIFIER 


DISCRIMI- 
NATOR 


ANTI- 
COINCIDENCE 
CIRCUIT 


Fig. 1. Block diagram of the liquid scintillation counting system. 


_ The sintillator cell and the photomultiplier tubes are placed inside 
a large anti-coincidence ring counter filled with 1 atm. propane and 
working in the proportional region. This is the same shielding 
system that is used for radiocarbon dating by Nydal4). The 
counters are surrounded by 20cm of iron, and the laboratory 
itself is situated 4m below ground level with a total of 1 m con- 
crete above. 

A scintillator solution consisting of 4g PPO and 0.1 g POPOP 
per liter toluene was used throughout these experiments. The 
multiplier tubes were kept in darkness for some weeks prior to and 
during the measurements. 
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§ 3. Radiocarbon counting characteristics. Following Hayes and 
Ott °), we have chosen the high voltage and the lower pulse height 
discriminating level of the counting system so as to maximize the 
value of E?/B, E being the radiocarbon counting efficiency and B 
the background counting rate. No upper discriminating level was 
used in this work, although some additional increase in E2/B 
probably could have been achieved by this means. 

The curves of fig. 2 illustrate the background counting rate and 
the radiocarbon counting efficiency as functions of the lower 
discriminating level at the optimum multiplier high voltage. The 
curves are measured with the anti-coincidence shielding system 
active. The optimum lower discriminating level is indicated. 


| 
| 
OPTIMUM. 
DISCRIMINATOR 
LEVEL 


V4 


BACKGROUND COUNT RATE,COUNTS/S 
& 
RADIOCARBON COUNTING EFFICIENCY; 


% 50 100 150 
DISCRIMINATOR LEVEL, ARBITRARY UNITS 


Fig. 2. Background count rate and radiocarbon counting efficiency as 
functions of the discriminator (lower) level at the optimum value of multiplier. 
BOE: 


The performance data of the scintillation counter operated at 
room temperature may be summarized as follows: 


Radiocarbon counting efficiency Sele 
Background, anti-coincidence shield passive 51 counts/min. 
Background, anti-coincidence shield active 16 counts/min. 


There are experimental indications that light interaction between 
the multipliers and radioactive contamination in the steel cylinder 
are responsible for the main part of the residual background. The 
radiocarbon counting efficiency was determined using a benzoic 
acid (carboxyl C-14) standard, supplied by The Radiochemical 
Centre, Amersham, with a stated accuracy of 5%. 
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§ 4. Activity of contemporary wood. The counter has been used to 
determine the radiocarbon contents of ethanol made from con- 
temporary wood. 7 ml of 96% ethanol was added to 93 ml of 
scintillator. Under the assumption that this would not change the 
radiocarbon counting efficiency (6), the specific activity was found, 
to be 

13.1 -- 0.6 disintegrations/min.gram carbon. 

This is in excellent agreement with recently reported values: 


Hayes andArnold’) 13.6 + 0.2 disint./min. gram carbon, 
Fergusson 8) 12.00 OL ee, Zs ee 
Nydal 9) 13.0;5=O9l on. ~ s ie 


The last value pertains to the National Bureau of Standards oxalic 
acid “contemporary carbon”’ standard. 
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A SIMPLE LOW BACKGROUND PHOTON 
COINCIDENCE DETECTOR FOR TRITIUM LIQUID 
SCINTILLATION COUNTING *) 


by R. S. SIGMOND and K. G. SCHJETNE **) 


Institutt for Teknisk Fysikk, Norges Tekniske Hogskole, Trondheim, Norway 


Summary 

The electronic system of this coincidence counter consists of two. 13-stage 
Venetian blind photomultipliers coupled to a simple DeBenedetti-Richings 
coincidence circuit, followed by a medium-gain amplifier, pulse selector and 
scaler. The main feature of the system is its simplicity, chiefly realized by 
utilizing the quasi-linear transfer properties of the take-the-lesser type 
coincidence circuit. ‘‘Balance-point operation’’ permits the use of multiplier 
high voltage supplies and pulse amplifiers not possessing the high stability 
usually required for scintillation counting. The counter has been operated 
at 26°C with a tritium counting efficiency of 16° at a background of 100 
c/min. The minimum detectable tritium concentration in water is 1.7 « 10-8 
curie/liter. However, the time needed for the counter to ‘‘quiet down”’ after 
a change of sample is prohibitively long, 2-3 days. Of the background, 
nearly 50% seems to be caused by light interaction between the photo- 
multipliers. Some suggestions for future improvements are given. 


§ 1. Introduction. The main object of this investigation was the 
development of an efficient low background photon detector for 
tritium liquid scintillation counting, complying with the accessory 
conditions of simple electronics and room temperature operation. 
The instrument was chiefly intended for the investigation of natural 
tritium in water. 

Several tritium liquid scintillation counters of fairly high ef- 
ficiency and low background have been described in the litera- 
ture 1-11), Both single photomultiplier tube spectrometers ”)?)4) and 
coincidence spectrometers have sucessfully been used for room 


*) Condensed from the N.T.H. Technical Report T.. 60-1. 
**) Now at JENER, Kjeller per Lillestrém, Norway. 
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temperature counting of tritiated water; the best performances 
being efficiencies of 10-20% at backgrounds of 100-200 counts/min 
for water sample sizes of about 1 ml. However, in this application 
the single tube spectrometers require multiplier tubes specially 
selected for low dark current noise, while the coincidence spectro- 
meters reported all feature a complicated and costly electronic 
system. 

This counter employs a coincidence photon detector and thus 
does not require doubly selected multiplier tubes. Nevertheless its 
electronic system is almost as simple as that of the single tube 
counter. This has been achieved by using very high gain multipliers 
coupled through cathode followers to a coincidence circuit whose 
quasi-linear transfer properties permit any pulse height selection 
to be done after the coincidence selection. The instrument does not 
have the linear energy response of the commercial scintillation 
spectrometers and the threshold energies are accordingly more 
cumbersome to adjust. However, this constitutes only a minor 
drawback when only one specific radioisotope is to be measured. 


CHANNEL 
ANALYZER 


Fig. 1. Block diagram of the liquid scintillation counter. 


§ 2. Photon detector electronic system. A block diagram of the 
liquid scintillation counter is given in fig. 1, while fig. 2 shows the 
circuit diagram of the photomultiplier-cathode follower-coincidence 
system together with a sketch of the scintillator cell and photo- 
multiplier housing. The photomultipliers chosen are a matched 


TRITIUM LIQUID SCINTILLATION COUNTER 95 


pair of EMI 9514A. This tube is of the Venetian blind type with 
13 dynodes and a S11 CsSbO photocathode, and can be used at 
gains exceeding 108. It has a transit time-spread of 16 nanoseconds, 
more than twice that of the best focused-dynode types, but this 
drawnback is offset by a somewhat lower dark current. Moreover, 
the following coincidence circuit can be comparatively slow, thus 
making it easier to obtain some degree of linearity in the coincident 
pulse amplitude transfer function. 
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Fig. 2. Diagram of the photomultiplier — cathode follower — coincidence 

circuits. Sketch shows scintillator cell, photomultiplier and cathode follower 
housing. 


The coincidence circuit is due to De Benedettiand Richings!%). 
It is of the take-the-lesser type, and its output pulse amplitudes 
thus contain a fair amount of information about the photon pulses 
emitted from the scintillator. The resolving time of the circuit, as 
defined by Bay13)14), is strongly dependent on the input pulse 
amplitudes, increasing from 20 ns for signals of amplitudes equal 
to twice the lower discriminating level, to 80 ns for the maximum 
signal amplitudes available from the cathode followers. Most of 
the uncorrelated noise pulses have amplitudes below twice the 
lower discriminating level. 

The output pulses from the coincidence circuit are quite slow, 
thanks to the pulse stretcher diode arrangement. Accordingly the 
following pulse amplifier, pulse height selector and scaler need 
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only have resolving times better than some tens of microseconds, 
as long as the expected counting rates are low. Furthermore, poor 
overload properties of the amplifier can be tolerated because the 
cathode followers preceding the coincidence circuit limit the 
signal amplitudes. In the present investigation the amplifier and 
pulse analyzer system had a resolving time of about one microsecond. 

Coincident and single signal response curves for the coincidence 
circuit, amplifier and pulse height selector system are given in 
fig. 3. The tritium channel and the maximum coincidence circuit 
input signal amplitude are indicated. It is seen that no single signals 
can enter the tritium channel. 
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Fig. 3. Coincident and single signal response curves for the coincidence 
circuit, amplifier and pulse analyser system. The tritium channel and the 
maximum available input pulse amplitude are indicated. The single signal 
response curve is shown with the vertical scale enlarged 10 x. 


At multiplier gains of 108 the tritium channel lower and upper 
levels are equivalent to 0.5 and 2.2 electrons leaving the photo- 
cathodes, respectively. Taking account of the statistical nature of 
the electron emission and multiplication processes and of the take- 
the-lesser operation of the coincidence circuit, this is just the main 
range of yields to be expected from tritium scintillations. 

The present coincidence system has worked very well and is not 
critical with respect to component values. The only parameter that 
has to be adjusted with some care is the coincidence diode forward 
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bias current. An increase in this current causes an extension of the 
range of linear coincidence response and a decrease in the single 
signal rejection ratio. 


§ 3. Lvitiwm counting *). The work reported here was mainly 
concerned with the development of the photon detector electronic 
system. Accordingly we have chosen a scintillator and an optical 
system believed to be fairly efficient, without bothering much about 
optical optimalization or ease of sample changing. 

The 50 ml cylindrical quartz scintillator cell has plane end 
windows optically coupled to the multiplier photocathodes with 
thin layers of immersion oil. The cylindrical part of the wall is 
coated with a cellulose-base TiO laquer. The scintillator solution 
used contains 45 ml of 50 g/l napthalene, 7 g/l PPO and 0.05 g/l 
POPOP in #-dioxane, mixed with 5 ml water sample, and bubbled 
with 99,99°% argon for 15 minutes in the cell 1) 45), The multiplier 
tubes had been kept in darkness with high voltage applied for some 
months prior to the measurements reported below, but had to be 
exposed to red darkroom light for some minutes during sample 
changing. The cell, photomultipliers and cathode followers are 
placed inside a steel-lined 5 cm lead shield. The operating temper- 
ature of the scintillator and photocathodes was about 26°C during 
the measurements. 

As the coincidence circuit has a very restricted range of linear 
response, the lower discriminating level as referred to the input 
must be considered as fixed. Thus the only important parameters 
that can be varied to make the tritium counting properties of the 
system optimum are the upper discriminating level and the multi- 
plier high voltage. 

The upper discriminating level has been set safely below the 
upper (overload) hump in the background pulse height spectrum 
(see fig. 4). The multiplier high voltage has been adjusted to 
maximize the tritium counting efficiency, thus making this im- 
portant quantity insensitive not only to small variations in the 
multiplier gain, but also (because of the quasi-linear transfer 
properties of the coincidence circuit) to drifts in the gain of the 


*) An earlier version of this counter has been used by Nygaard 2%) for radiocarbon 


measurements. 
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pulse amplifier. This balance-point-operation 16) will of course 
only apply to one specific isotope. 

Fig. 4 shows the tritium and background pulse height spectra as 
seen by the pulse selector under the above conditions. Finally, 
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Fig. 4. Background and net tritium pulse height spectra. The tritium 
channel is indicated. 
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——= DAYS AFTER FILLING OF SCINTILLATOR CELL 


Fig. 5. Mean total count rate in tritium channel for eight 5 ml tritiated water 

samples, activity 104 disint./min, as function of time after filling of the 

scintillator cell. The shaded area represents the total spread of the eight 
count rate curves. 
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fig..5 and fig. 6 show mean tritium sample and background count 
rates in the tritium channel for several samples as functions of time 
after filling of the scintillator cell. The shaded areas represent the 
total spread of the count rate curves. It is seen that some samples 
needed as much as three days to establish constant count rates. 
We have not investigated this effect further, but it might be due 
to the high operating temperature of the scintillator. However, 
after three days all the measurements on the 5 background and the 
8 tritium samples agreed with the following counter performance 
data within the stated standard deviations: 
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—»> DAYS AFTER FILLING OF SCINTILLATOR CELL 


——= BACKGROUND COUNT RATE, PULSE /MIN 


Fig. 6. Mean total count rate in tritium channel for five 5 ml distilled water 
samples (background) as function of time after filling of the scintillator cell. 
The shaded area represents the total spread of the five count rate curves. 


Background count rate, tritium channel, 26°C: 100 c/min + 7%, 

Tritium counting efficiency, tritium channel: (15.7 + 1,4)%. 
The 5% uncertainty in our Harwell tritium standard has not been 
included in the efficiency standard deviation given above. 

Following standard radiocarbon dating practice we may define 
the minimum quantity of tritium that can be detected as that 
giving a standard deviation of 25%. The sample and background 
counting times need not be specified in our case, the limiting 
factors obviously being the reproducibility of the background, 
and, to a lesser degree, that of the tritium counting efficiency. We 
find: 

Minimum detectable amount of tritium: 8.6 x 10-14 Curie, 
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Minimum detectable tritium concentration in water: 1.7 x 1078 
Curie/liter. 

Some measurements were made of the effects of the optical 
coupling oil and the TiO diffuse reflector on the tritium counting 
efficiency: 

Air coupling, no reflector: 2%; 
Oil coupling, no reflector: Ay fo 
Air coupling, TiOg reflector: 13%, 
Oil coupling, TiO reflector: 15.7%. 


The importance of good photon collection efficiency is evident. 

A rough experimental analysis of the background gave the 
following approximate results: 

Accidental coincidences between uncorrelated multiplher dark 
current noise pulses: 40 c/min (30°C), 20 c/min (26°C), 15 c/min. 
(24°C) 10'cjaman (2126); 

Coincidences due to light interaction between the two multipliers: 
_~> 45 c/min. 


§ 4. Conclusions. The coincidence counter for liquid scintillation 
counting of tritiated water described here has some notable ad- 
vantages: Room temperature operation, high concentration 
sensitivity and simple electronic system. Balance point operation 
with respect to variations in the gain of the electronic system 
permits the use of photomultiplier high voltage supplies and pulse 
amplifiers not possessing the high stability usually required for 
scintillation counting. 

But, in the present stage of development the counter also has 
some serious drawbacks, the most notable apperantly being 
associated with the scintillator. The time usually needed for the 
counter to “quiet down”’ after a change of sample is prohibitively 
long, 2-3 days, and even after this period the reproducibilities of 
the tritium counting efficiency and the background are relatively 
poor. . 

As for the future work to be done on this counter, there is a most 
obvious need for an investigation into the erratic behaviour of the 
scintillator. Furthermore, light interaction must be minimized 
by mounting the multipliers at right angles to each other. Stabiliza- 
tion of the operating temperature by means of tap-water cooling is 
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clearly desirable. The cathode followers preceding the coincidence 
circuit may be omitted. 
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MEASUREMENTS ON THE MAGNETIC ANISOTROPY 
IN EVAPORATED IRON FILMS 


by R. VRAMBOUT and L. DE GREVE 


Instituut voor Lage Temparaturen en Technische Fysica, Leuven, Belgium 


Summary 


We have examined the problem of the magnetic anisotropy in evaporated 
iron films by use of the magneto-optic Kerr effect method. As a result we 
found that there exists a weak anisotropy in some iron films coming from 
mechanical tensions. 


§ 1. Introduction. Since a few years we have been especially 
interested in the magnetic structure of evaporated iron films 1)?). 
At present we report on the result of our latest measurements on 
magnetic anisotropy, obtained with the use of the magneto-optic 
Kerr effect method. This method does not only allow to study the 
hysteresis phenomenon of the whole film, but also at any particular 
point of the film. Indeed, the illuminated surface of the film can 
easily be reduced to a few tenths of a millimeter. 

It may now safely be assumed that thin ferromagnetic films 
are composed of a single Weiss domain, especially after being sub- 
jected to a magnetic field. Before magnetization, Bloch walls may 
in some cases be found 1), but as a rule they disappear completely 
when an external field is applied. 

We have tried to determine the different factors that influence 
the magnetic properties of these one-domain iron layers. When 
measuring as a function of thickness and especially when evaporat- 
ing in a rather poor vacuum, we experienced earlier that it was 
practically impossible to obtain reproducible results. So we decided 
to restrict our investigation to layers of practically the same 
thickness (between 730 and 1120 A). According to Kittel’s theory 3) 
and our own previous measurements these films should indeed be 
composed of only one domain. 
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§ 2. Experimental method. Our experimental method has been 
briefly described in an earlier paper 2). In the Kerr effect a polarized 
lightbeam is reflected from the layer that is to be studied. This 
causes a small rotation of the plane of polarization. Our actual 
experimental device differs only little from the previously described 
apparatus. The CdS photoconductive cell, that was then used in 
order to measure the changes in light intensity after reflection, 
was of course not sensitive enough to measure differences in 
coercive fields with a sufficient accuracy. So it was replaced by 
a photomultiplier, in this way increasing the sensitivity of the 
method by several orders. This allowed us to use a weaker source 
of light. One of the disadvantages of a high pressure mercury arc 
lies in the difficulty of getting it sufficiently stabilised. So the 
mercury arc was replaced by a much weaker incandescent lamp of 
the type commonly used in microscopes. 

The iron film can be rotated in its own plane between two Helm- 
holtz coils in such a way that the magnetic field always remains 
parallel to the plane of the film. Hysteresis loops are then measured 
in different directions of the film. The angle between the long axis 
of the rectangular glass support and the external magnetic field 
varies between 0° and 180°. 

The films were evaporated in a vacuum of the order of 10-5 to 
10-6 mm Hg on polished rectangular glass plates. The dimensions 
of the layers were of the order of 3 to 4 mm in either direction. The 
thickness of the layers was measured by means of the optical 
multiple interference method of Tolansky. 


§ 3. Experimental results. The determination of hysteresis loops 
as a function of the position of the layer is a very reliable method 
for detecting the existence of any anisotropy axis. Indeed, any 
change of magnetization in such films can only be attributed to a 
rotation of the magnetic moment in the whole film and not to the 
displacement of Bloch walls. Hysteresis loops which are measured 
parallel to an anisotropy axis should be rectangular, while loops 
measured perpendicular to this axis can be expected to be almost 
flattened out. From our measurements it appeared that often a 
relatively weak anisotropy occurs. Nevertheless in some layers no 
anisotropy can be detected. Mostly, however, the shape of the 
hysteresis loop and also the coercive force changes as a function 
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of the external magnetic field direction. Fig. la and 16 reproduce 
two hysteresis loops measured respectively parallel and perpen- 
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Fig. 1. Hysteresis loops measured 
a) parallel b) perpendicular to the long axis of the glass support. 


dicular to the long axis of the glass support for film 11. The different 
shapes of the curves can clearly be noticed. Fig. 2 gives the variation 
of coercive force as a function of the angle. There remains no 


H(8)) 
%0 


x 0,383 oersted ——> 
3 
=} 


a) 
Ss 


O° 225° 45° 675° 90° N2s° 135° {f5%S° io° 


—_— 


Fig. 2. Variation of coercive force as a function of the angle. 


doubt that there exists a weak anisotropy axis lying closely to 
the long axis of the glass support. 

The direction of the anisotropy axis corresponds to the direction 
of the maximum coercive force. As a matter of fact the anisotropy 
axis if existing is in all our layers either parallel to the long side of 
the glass support or very close to this direction. Table I contains 
the thickness D of the films, the angle @ between the anisotropy 
axis and the long axis of the support, the magnitude H, of the 
coercive force measured parallel to the anisotropy axis and its 
maximum variation A H,, i.e. the difference between the values 
of the coercive force measured either parallel of perpendicular to 
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No. Film 


AH, (oersted) 


OnN anf WN 


TABLE I 
D(A) 0° H, (oersted) 
860 alee 58 
830 na. 49 
890 n.a. 14 
890 n.a. 25 
730 Of ced 
730 O° 34.4 
730 Pee 40.2 
760 n.a. 63.2 
760 OS 43.3 
760 n.a. 42.8 
745 10° 53.6 
1120 4° 73.9 


n.a.: no anisotropy could be detected. 


5 


the anisotropy axis. We must notice that the anisotropy field AH, of 
all films is small compared to the intrinsic coercive force (H — AH) 


of the film. 


Hysteresis loops have also been measured at particular points 
in the layer by cutting down the illuminated area to a few tenths 
ofa millimeter. Fig. 3a en 3d illustrate the hysteresis loops of film 11, 


a 


10,383 oersted 


light intensity 


b 


Fig. 3. Hysteresis loops for a small area of the film measured 
a) parallel 6) perpendicular to the long axis of the glass support. 


x0,383 oersted 


but this time measured at a point of the film. These hysteresis loops 
are somewhat more rectangular than hysteresis loops of the film 
as a whole. The differences AH, between maximum and minimum 
value of the coercive force slightly increases (6.9 instead of 6.1). 
In other films or at other points A H; can also be smaller. This 


106 MAGNETIC ANISOTROPY IN IRON FILMS 


means that these layers are slightly inhomogeneous, and that the 
values measured over the whole film can only be regarded as 
average values. 

The influence of the shape of the film on its magnetic properties 
was also examined. During the same operation we evaporated 
circular, square and rectangular films on equal rectangular glass 
supports. No essential difference in their magnetical properties 
could be found. This seems to indicate that any anisotropy in these 
layers should merely be attributed to mechanical tensions, especially 
if no crystalline anisotropy can reasonably be expected in such 
complicated polycrystalline structures. 
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NEAR-ZONE BACK-SCATTERING FROM LARGE 
SPHERES *) 


by V. H. WESTON 


The University of Michigan, Ann Arbor, Michigan, U.S.A. 


Summary 


For an incident electromagnetic plane wave, the near-zone behaviour of 
the backscattered field produced by a perfectly conducting sphere is in- 
vestigated for small wavelengths. The backscattered cross-section becomes 
appreciably different when the receiver approaches to within a distance of 
several radii from the center of the sphere, and in fact becomes the cross- 
section of a flate plate for the receiver very near the sphere. 


§ 1. Introduction. An abundance of work has been done on the 
various problems connected with electromagnetic scattering from 
a sphere. Asymptotic results for small wavelengths have been given 
for forward, back and total cross sections. However, apart from 
results giving the currents on the sphere or the local behaviour 
of the fields in the transition region of the geometrical shadow 
boundary, the results hold generally only for the far field. 

‘However, there arise many cases wherein the near field is impor- 
tant. One can see this in an example ofa satellite orbiting the earth, 
or a planet at a distance of a radius from the surface. The satellite 
certainly is in the near zone of any signal bounced off this planet. 

Here the problem that will be considered is that of the near field 
behaviour of the back-scattered field produced by a high frequency 
electromagnetic plane wave incident on a perfectly conducting 
sphere. It will be shown that the back-scattered cross-section differs 
considerably on approaching close to the scattering body. 

In § 2a modified Watson transformation will be performed on the 
exact Mie series answer for the back-scattered field. The behaviour 


*) The research reported in this paper was supported under contract AF-30 (602)-1853 
with the Rome Air Development Center. 
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of the portion of the scattered field due to the diffracted waves 
will be investigated for the near zone. In §3 the asymptotic ex- 
pression will be deduced for the portion of the back scattered field 
corresponding to the reflected wave for the far field. In § 4 the re- 
flected field portion will be obtained for the near field using the 
Kline-Luneberg expansion. 


§ 2. Behaviour of the diffracted field contribution. For a plane 
harmonic wave of the form 


Et oe te elke tet (1) 


incident on an large perfectly conducting sphere of radius a, the 
back scattered field is given by (ref. 1), p. 564) 


co 


Es = —ig ett  (—i)" (n + 8) (an + ibn), (2) 
where 
bah hn (RR) Jn(ka) 
| Pee eae Ss 
and 
ee [ARAn® (RkR)]’ [kajn(Ra)]’ (4) 


kR[ Rahn“ (ka) |’ 


Henceforth, the time factor e~'” will be dropped. Using the series 
expansion 


e-tkR _. 9 E (i) (n + 4) jn(RR), (5) 


the following is obtained: 


E*(k) = (I, 4 1f2) + a eikR—2ika __ at eWKR __ gtk: (6) 
where 
fr = B (—1)"(n + 3) [in(RR) — an())- (7) 
and 
os [ [RRjn(kR)]’ | 


Iz = ¥ (—i)™ (n + 2) aR n(h) (8) 


n=0 
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The two series may be summed by contour integration giving 


: a ‘ ) a ie : E (RR 
Ace 2kR) 4 COS vz Poa a 


C 


oe ] du (9a) 
H,) (ka) 
and 
Ip = ee: a eae ” VER Hy kR)]’ — 
: 2/ 4kR J cosvx its Sees) 


C 


[VAR Hy (kR)]’ 
[V ka Hy) (ka)]’ 


[Vka H,2 (ayy clonme F (95) 


where C is the contour surrounding the poles v = 4 
shown in fig. 1. 


3, 3, ... as 


’ 


v plane 


Fig. 1. Contour C in the v-plane, surrounding the poles v = §, 3, 8, .... 


It can be shown that 


sama (10) 


(o. Dez Di 


where Dz is the contour surrounding the poles of the integrand 
(i.e., for Ty the zeros of H,“)(ka), and for Ig the zeros of 


(Vka Hy) (ka)]’) 


and D, isa line integral extending through the origin of the v plane. 
The line integral can be replaced by the integral containing the 


110 V. H. WESTON 


even part with respect to v of the integrand, and this in turn can 
be replaced by twice the integral extending from the origin to 


infinity, giving 
oo exp(—7f) 


ie eae ae | F sin vx 
2kR COS U7 
[Hy (RR) Hy (ka) — Hy (RR) Hy) (ka)] q 
2 ce?" (ka) 
45 er ie oe A, (RR) Hy, (Ra) 


p> 1291" 
2kR Pa 1+ e* d H,,,0 (ka) 
dd; 


v 


(11) 


where @, are the zeros of 
beR 1 
H, (ka) = 0, 5) = ka + ($ha)* et" G + 0 (aa) 


gy are the zeros of the Airy integral, and 0 < B < 32 
First consider the residue portion of (11). The dominant part is 


the term 
etvint A (RR). (12) 


Since d; lie above the real axis, e’?'* is a decaying exponential. The 
decaying exponent is of the order of 9:(ka)*. In the far field, ex- 
pression (12) behaves like 


etaint / 2 eikR-Moin—hiz 
: akR ; 


and in the region where RR > 1 and (ka)2/R ~ O(1), expression(12) 


y= exp[?djn — din + i(kR + ka2/2R) + 1(a/R)(kka)* e** Gi], 
and in the very near field as R — a, it becomes the order of 
eftai7 


A similar analysis holds for the residues obtained by transforming 
the contour integral representation for the series Ig. The residues 
corresponding to the diffracted field are obtained from the zeros of 


[V (ka) Hy (ka)]’, 


which are related to the zeros of the derivative of the Airy integral. 
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Essentially for very large spheres, the residue terms are more 
important in the far field than near field in the direction of back- 
scattering. Physically this is what is expected. To reach a point in 
the near zone, the diffracted waves must creep a greater distance 
around the sphere, and since their attenuation depends upon 
distance travelled as a surface wave, the diffracted waves repre- 
sented by the residue terms must be attenuated more in the near 
zone. 


§2. Behaviour of the reflected field contribution. The integral 
portion of expression (11) will now be considered. This can be 
calculated in the same manner as is done by Scott 2). First it 
should be noted that on letting 6 = 41, the integral 


Pe ci) : 
(se) 2 | y————_-#™ H,® (BR) do, 
0 
is identically 
4 eke. (13) 


To obtain an asymptotic expression for the second integral in 
expression (11), the asymptotic expression for large ka, 


an Hv) (RR) Hy) (Ra) ( 2 ‘8 


Ge hi = | ——_ 
H,® (ka) akR 


7 ( l l ) | sl 
- exp E — 2ika — wW ie aR raat 
is substituted in the integral, giving 
coexp(—if) 


1 1 ) 
; 2 fas oy gael See eee ng Ee 
v tan vx expi| RR — 2ka — v ( P PR 


1] 


2(kR) 
0 


—ivm 
and this, in turn, after expanding tan vz in a power series of e 


gives asymptotically 


= ae} exp i(kR — 2ka) [: +0 ai (14) 


The integral corresponding to Jz may be evaluated in the same 
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manner to give 


Le—ik wa aes 
de—*R (_; + 1/(RR)] 4 2(2R — a) E ; o(=_) |e ars) 


Hence, the reflected field portion of the back-scattered field is 


given by 
a ] 
_ 3. at(kR—2ka) (e230 ( I. 16 
is (2R — a) [ core Oy 


For more details concerning the analysis, the reader is refered to 
Scott’s paper. However, this technique is very difficult to use in 
obtaining higher order terms in ka for the near field. The near field 
representation for the higher order terms in ka will be obtained by 
the method given below. 


§ 3. Near zone reflected field. To obtain the near zone component 
for the reflected field portion, the Kline-Luneberg expansion will 
be used. 

When a plane harmonic wave is incident on a perfectly conducting 
body, the portion of the scattered field due to reflection only can 
be represented in the form 


Es = e'**(Eo + E,A + Eod? + ...). (17) 
If the incident radiation is being propagated in the direction of the 
z-axis and is polarized in the direction of the x-axis, and if the 
scattering body is a body of revolution with the z-axis as symmetry 


axis, then by Schensted 8) the above field components satisfy 
the relations 


} hyih 
= (es role 2 | 
n Iho te re A V Bass ds (18) 
f 


“ 
Qn 
where z = /(p) defines the surface of the body, and p = Vx2 + y?, 

The curvilinear coordinate system (p, s, $) is described as follows: 
Through each point in space (exterior to the body) there passes a 
ray which has been reflected from the body: 

p is the distance from the z-axis at which the ray hits the body ; 


and 


Eno = ——(V:En i(APa+ sy 2 AO, (19) 
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s is the sum of /(p), and the distance along the ray from the body 
to the point; 

¢ is the angle between the «z-plane and the plane formed by the 
ray and the z-axis. 
The metric coefficients are defined as follows: 


dl? = hy? dp2 + ds? + ho? d¢?; 


and 
Blass) 2f'(s — f) 
hy il iS : he ad + : —-, 20 
(41 He Parr ee 
In addition, Schensted 3) showed that 
Ey = ante (— cos dp + sin # 4). (21) 
jhe 


The problem, at present, is to evaluate the coefficients E, in the 

backscattered direction for the case when the body of revolution 

is the sphere. ; 
For a sphere of radius a, the equation of the illuminated face is 


i(p) = a — Va? — p?, (22) 


and the back-scattered direction represented in the coordinate 
system (p, s, ¢) is given by p = O. If the fact is used that the back- 
scattered field will be polarized in the same direction as the incident 
field, one obtains 


(En-iy), 0 = 0, (En-iz),—o = 0. (23) 
Thus one need only evaluate 
(En-iz)o—o: (24) 
The term Eo-i, can be expanded in the form 
Eo-iz = Ao + p2(Bo cos 26 + Co) + p4 (Do cos 26+ Eo) +..., (25) 


and it can be shown that 
En-iz = An+ p?(Bncos 2¢ + Cn) + p4(Dncos 26 + En) +... (26) 
Hence, for back-scattering 

Ad, Ain 


will be required, but for present purposes Ao, Ai, and Ag will be 
computed. 
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Now it can be shown that . 
E29 iz ~ O(p?), (27) 


and 
E,°-iz ~ O(p4). (28) 
Thus, one obtains 
1 hyhe 


‘oe _[ VB pace to as + 016, (29) 
47 p P 
0 


and 


; hyhg Iyh 
Ey -iy = i / : {V ck is V2(Eo-iz) ds + O(p4). (30) 
7 ms p 


In order to compute (Eg-iz), 9, the term [V?(E1-iz)],9 is re- 
quired. 
It can be shown that if 


U= (Ay + p*Cy + pin + .--) 4 COS 20( Bap" 4 Dap =.=, alow) 
where An, Cn, Bn, ... are functions of s only; then 
V2u = (An* + p?Cn*+ptEn* +...)+c0s 26(Bn*p?+Dn*pt+..., (32) 


where 


An ( 4 Vee 


a 2 
n= — | —————]}, (33 
2 TAPS +a Os BEE ( ) (88) 


2s-+ a) 
and C,* depends upon An, Cn, En only. All the starred coefficients 
are functions of s only. 

Thus in order to compute [V?(E1-iz)], 9, only the coefficients Ay 
and Cy are required and, similarly, in order to compute A; and C; 
the quantities Ao, Co, and Eo are required. These in turn are found 
from the expression 


Eo:ig = en (— cos dp + sin ¢). -iz 
Wg 


$id a eel if) Aube lope 
eed ee ee 
= (Ao + Cop? + Eop* +...) + cos 26(Bop? + ...). (34) 
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Computing the coefficients, the following is obtained: 


Erie eee a! 
A 1 l 2a 2 
Br tao =| Ta 
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f —] 
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Now for back-scattering, s is given by the relation s = R — a, 
where R is the distance of the observation point from the center 
of the sphere. Hence, the reflected field is given by 


: a a ay a 
Boe Ro a [ SAR a 
a3 
gaat al (38) 
where 
Pia sina (39) 
(QR =a)?” 
and 
yas a(R — a)(2R? — 4Ra + 3a?) (40) 


(2k — a)4 


§ 4. Comments. The back-scattered cross-section is defined by 


o = lim 427? |E$|2/|E*/2. (41) 
+00 
However, this is usually defined for the far-field where the distance 
rv is the distance from the scattering object to the receiver. However, 
in the near zone the distance 7 is ambiguous since the body no longer 
looks like a point source. For near-zone back-scattering from the 
sphere, it is best to take 7 as the distance from the specular point 
to the receiver. Hence, if R is the distance from the centre of the 
sphere of radius a, 7 is given by ry = R —a. In the far field 7 
asymptotically approaches R; hence, in the cross-section expression, 
y may be replaced by R. 
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For the near field the back-scattered cross-section may be obtained 


using (38), giving 


o = 02 ¢ as a I +0 (az) | (42) 


Hence, it is seen that the cross-section decreases as the receiver 
approaches the sphere. These results hold, provided that the 
wavelength is sufficiently large, such that the diffracted field is 
negligible in the far zone. It was shown in § | that the contribution 
of the diffracted field to back-scattering is even less in the near 


field. 


Note added in proef: 
It is interesting to note that 
G = Anrt: ar ea. 
C= ara’, Y= "4; 
C= TA. ToS a. 
For the first case, the backscattered cross-section is that of a flat 


plate which is expected, and the third case is the usual far field radar 
cross-section of a sphere. 
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PIELDS IN-SQUARE HELMHOLTZ COILS 
by R. D. STRATTAN and F. J. YOUNG 


Carnegie Institute of Technology, Pittsburgh, Pennsylvania, U.S.A. 


Summary 

The axial component of the magnetic field in square Helmholtz coils 
of rectangular cross-section is investigated. Expressions for the axial field 
are derived and simplified enough so that any specific case can be easily 
calculated with a small internally programmed digital computer. Plots 
of the field at the centre of the coils as a function of winding thickness with 
coil spacing as a parameter are presented. Curves of the variation of the field 
with distance from the centre are given. Spherical volumes in which the 
axial magnetic field varies 1% and 2% are computed for various con- 
figurations. The dimensions of the coils which hold the field variations to 
1% and 2% over the maximum spherical volume are given. 


§ 1. Introduction. The square Helmholtz pair shown in fig. 2 
can be used to produce uniform magnetic fields in the vicinity of 


/ “ 


Fig. 1. Square current filaments. 


at 


the origin of the coordinates. In applications such as the rictometer 
where more than one Helmholtz pair is needed square Helmholtz 


ei lf — 
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coils are more compact and consume less power than ordinary 
Helmholtz pairs. In the case of 6 = 0.80 and «= 0.15 for the 
inside pair of a two phase-arrangement consisting of two Helmholtz 
pairs a circular set of coils require about 43% more copper to 
produce the same magnetic field at the same current density than 
a set comprising square coils. The copper losses are correspondingly 
ereater for circular than square coils in arrangements where more 
than two Helmholtz pairs are required. 


§ 2. Calculation of the magnetic field. The field due to a coil 
system similar to the one shown in fig. 2 can be found at any point 


Fig. 2. Square coils filling a cubical. 


in space from Ampere’s law and integration over the winding 
cross-section. The final integral cannot be performed analytically 
but is easily evaluated with the aid of a digital computer and 
Simpson’s rule. The only case calculated here is for the overall 
axial dimension equal to the transverse dimensions. However this 
method is applicable to any rectangular prism shape. 

The parameters of this coil system are a, the dimensionless 
spacing between the coil sections, and f, the dimensionless measure 
of the winding thickness. The gap of 2«L is placed between the coil 
sections to give a more uniform field and to provide better access 
to the interior of the coils. 

Only the axial field component is calculated, although the radial 
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components can be calculated by the same method. First the axial 
field component due to the two current filaments shown in fig. 1 
is found from Ampere’s law. Then using this as the field due to the 
infinitesimal winding areas dd by dé, the total field is found by 
integrating over d from «L to L and over ?¢ from pL to L. Rationalized: 
MKS units are used in this paper. 


5+- 


4nHAL AT ORIGIN 
Np 


ri) 6 Us 8 : 1.0 


DIMENSIONLESS WINDING 
THICKNESS B 


Fig. 3. The magnetic field at the origin as a function of coil spacing and 
winding thickness. 


The field dH at (x, y, z) due to the current elements dl; and dl 
in fig. 1 is given by 
dH = I(dk, x a,,)/4ary2 4+ I(dlz X a,,)/4a7r2?, 
as a,(x — x1) + a,(y + d) + a,(z — 21) 
re [we —m)® + (y+ a)? + @ — 2)! 
a a,(x — x2) + ay(y — d) + az(z — 22) 
re [(e — xe)? + (y — a)? + (z — 22)" 
X41, %2, 2, and z2 = +f in order to define the current filament. 
Taking only the y component of the field, 


AnH -{{ («+ #) des 

Io JS ((e +82 + (y + 4)? + (2 — 21)?]? 
a (x + t) dzs 
le +92 + yy — a? + @— 2) 


a 
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(¢ — x) day 
[(e — 8)? + (y+ dy? +e — a) 
(¢ — x) dze 
(4 — 1)? 4 (y= a)? agi 
(z + ) dxy 
" [(@ — ma)? + y+ a? + (@ + HI 
(z + t) dx 
[(% — x2)? + (y — 4)? + (2 + 27]? 
(¢ — z) dx 
[( =a)? ae pie Oy? alee ae 
(¢ — 2) de ! 
* Te — m2 +o —4?+ 0a 


This integral can be evaluated using the following identity: 
t 


{ dx “i 1 | (¢ + a) 
(ere Tae +e) iat ort e+ oh 
= 

(¢ — a) | 
[a — a? + be + ort J 
If we define a quantity C;,;,4(x/t, y/t, z/t, d/t) as follows; 

1 

MGS ltl? eae t= 1 alee 


a 


Cin = 


1 
[1 + (—1)F2/t}? + [d/é + earn 
[1 + (—1éx/altt + (—1)%2/E] 
[(1 + (—1)t/é)? + @/é + (—1)sy/t)? + (1 + (1) e/#) 2] 
then the field due to the current filaments of fig. 1 is 
4nH i hy 
Dy) D5 y (1/t) Ci,j,n (x/t, y/t, 2/t, d/t). 


= Wiha eal il Aa | 


The field due to the coil configuration shown in fig. 2 is then 


Lit 
4nH 


1 
Yor Rages 
TL =z = faq) [acai Ci,g,n(x/t, y/t, z/t, d/t). 
B 


aL /t 
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With the substitutions 
[(1 + (—1)t&/é)2 + (1 + (—1)*2/é)2]}# tan uw = d/t + (—1)4y/t 


and sin « = w the integration over d/t reduces to the form 


rdw 
| p2 + w? = (1/p) arctan (w/f). 


4H ALONG Y-AXIS IN % 


Es 


—— SQUARE COIL 
—-- CIRCULAR COIL 
o<=15 , B=.80 
(b) 


OH ALONG Y-AXIS IN % 


OH ALONG Y-AXIS IN% 


Fig. 4. Variation of axial field uniformity along the y-axis; 
(4) 8. = 0275, (b) B = 0.80, (2) pi—0180: 


The result is 


Veo d(d/t) = 5 3 (—1)™ arctan 7%, (x/L, y/L, 2/L, t/L) 


aL/t 


m=1 n=1 
where for m= 1 and n = 1 ¢$f7% is defined as 
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[e/L + (—1)*2/L)la + (—1y/L] 
[E+ (—1)x/LVEIL + i 
4 (= 1) JL)? Ey (1) Rei)? lor Ny 
For m = 2, the expression is the same except that a is replaced by 
one; and for n = 2, (—1)#x/L and (—1)’z/L are interchanged. The 


total field in the axial direction then becomes 
1 


teen a) eke 
=— = Tia » 2 l( 1ym [arctan CMR ENE SIE cs te, 
ie j=1 j=1 k=1 m=1 n=1 i/L) a(t/L). 


ZA H ALONG X-AXIS IN % 


ZH ALONG X-AXIS IN % 


H ALONG X-AXIS IN % 


4a 
b 


Fig. 5. Variation of axial field uniformity along the x-axis; 
(a) B =0.75,. (b) B = 0.80, (c) B = 0.85. 


This integral was evaluated numerically on the C.I.T. IBM 650 
digital computer using arctangent and Simpon’s rule subroutines. 
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§ 3. Results and conclusions. The field at the centre of the coil 
system as a function of « and f is shown in fig. 3. For the ranges of 
« and f shown the field at the centre is given very nearly by 
H = 0.716 (1 — a)(1 — f£) JL amp-turns/meter. 

The percent deviation of the field from its value at the origin is 
given at points along the y and x axes by figs. 4 and 5 respectively 
for several values of dimensionless winding thickness. 


VARIATION IN H 


2.= ALONG X-AXIS 
X= ALONG Y-AXIS 


MAXIMUM DISTANCE FROM ORIGIN 
FOR I% 


VARIATION IN H 


FOR 2% 


A= ALONG X- AXIS 
X= ALONG Y-AXIS 


MAXIMUM DISTANCE FROM ORIGIN 


20 i 
DIMENSIONLESS SPACING o€ 


Fig. 6. The maximum distance from the origin for (a) 1%; or (b) 2% field 

variation. Intersection of equal f curves yield values of « and f for spherical 

region of uniformity. Maximum spherical volumes occur for a = 0.185, 

B = 0.80 with radius = 0.38 L at 1% uniformity and WO:45ulb wate 26 
uniformity. 


The maximum distance from the origin for 1% and 2% variation 
in the axial field as a function of « for a given f is shown in fig. 6. 
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The abscissae of the intersections of equal f curves are the «’s for 
maximum radius of a spherical volume of 1% or 2% uniformity. 
The maximum spherical volume occurs for 6 = 0.80, « = 0.185 
with a radius for 1% variation of 0.38 L. 

The percentage variation of the axial field of a solenoid of length 
2 R, outer radius R, inner radius 0.80 Rk with the windings omitted 
for 0.15 RF on either side of the midplane is shown in fig. 4(b) for 
comparison. The field at the centre of the solenoid is H = 0.115/R 
amp-turns/meter, while that of the cubical coil is H = 0.122]L 
amp-turns/meter for « = 0.15, B = 0.80. 


Received 17th October, 1960. 
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ASPECTS OF ELECTRICAL BREAKDOWN OF 
LIQUID INSULATING MATERIAL II 


by J. A. KOK and C. E. G. M. M. VAN VROONHOVEN 


Philips Research Laboratories N.V. Philips’ Gloeilampenfabrieken, 
Eindhoven-Netherlands 


Summary 


Part I of this investigation is continued by considerations on the ad- 
sorption of soaps and resins to the impurities which float in the in- 
sulating liquid. It was shown that soaps may stabilize the suspension, and 
resins may do likewise, but if applied both in a certain proportion immediate 
flocculation and lowering of the breakdown-strength may follow. The action 
of colloid chemical stabilizers or inhibitors originating from research on 
dyes has been tentatively explained. 


$1. §1. Introduction. In part I (ref. 1)) a discussion was given 
of the relation between the breakdown-strength of liquid insulators 
and the diameter of the insoluble impurities. If the radii 7 of the 
particles are equal and if their dielectric constant is infinite, then 
the polarizability may be assumed to be proportional to 73. If there 
happens to be a local stress concentration of thrice the average 
field strength, such as occurs on top of a hemispherical hump on 
one of the electrodes, then the breakdown strength Eo, occurring 
after a long duration of time, may be shown to depend only on 7, 
the absolute temperature 7, Boltzmann’s constant k and the 
permittivity of the dielectric eo: 
E973 E 92 te tRT. (1) 
With T = 300°K and «9 =1 this formula gives for Eo = 
= 25kV/mm:r = 25 A; Ey = 100 kV/mm:7 = 9.5 A. The relation 
was shown to be valid for suspensions of particles having a diameter 


of about 1 micron (10000 A). 
The size of the particles in aqueous suspension may be treated 
by means of considerations based on calculations carried out by 


— 125 — 
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Overbeek, Verwey and Hamaker 2). By the attractive London- 
van der Waals forces the particles tend to unite. In an ionogenic 
medium such as water the particles will gather a repulsive ion 
atmosphere. Their repulsion is proportional to the particle radius, 
whereas the attraction by Van der Waals forces does not keep 
up with it. In aqueous media a suspension will thus be stabilized 
after the particles have obtained a radius of 500 to 3000 A. 

In pure hydrocarbon oil the stabilizing ion atmospheres of the 
particles are very diluted and may be stripped off by the electro- 
phoretic motion by collisions against ions. However, stabilization 
of a suspension in hydrocarbon oil may nevertheless take place if 
surface-active substances such as soaps are present in the oil. A 
well-known application of soaps in aqueous suspensions consists in 
keeping soil particles dispersed in the process of laundering, thereby 
preventing flocculation of the particles and graying of the laundry. 

Constructive criticism of the views expressed in the first article 
(I) was brought forward by H. Koelmans, who pointed at the 
difference in the mechanisms of stabilization in aqueous suspension 
and suspensions in several hydrocarbon oils. 


§ 2. Colloid particles in oil, stabilized by soluble soaps or oxidation 
vesins. Soaps in mineral oil are considered to lower its insulating 
value. It will be shown here that addition of a small amount of 
soaps may improve the breakdown strength of the oil by stabilization 
of the colloid impurities it contains. However, if oxidation resins 
are present in the oil, an increasing soap content may be the cause 
of flocculation and a decreasing breakdown strength. 

A simple fatty acid may be represented by RCOOH, in which R 
is an aliphatic chain. The corresponding soaps formed may be RCOO 
Na or (RCOO)2 Ca. The heavier soap molecules with 14 to 22 
carbon atoms may have a length of 20 to 35 A. Micellar structures 
around impurities in oil containing soap molecules, which tend to 
be adsorbed to the impurities, may easily obtain a diameter of 50 
to 100 A, which is the size we need for our theory of electrical 
breakdown. Physical adsorption (by mirror image forces or Van 
der Waals forces, etc.) or chemisorption may play a part in this 
surface activity. If the energy of adsorption is very large, chemi- 
sorption plays a part, whereas physical adsorption liberates a lesser 
amount of energy (see ref. 3), 
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As interfaces are the seat of electrical charges or electrical double 
sheaths, polarizable molecules or charged ions and particles may 
become adsorbed at the interface of different phases (see fig. 1). 
Fig. 2 shows several examples of soap molecules adsorbed to water 
droplets in oil, oil droplets in water etc. Sometimes a double cover 
of soap molecules occurs at a water droplet in oil. If such a double 
cover of soap molecules comes into existence, the diameter of the 


Water-Soap Soap-Oil 


WLS 


Water 


Water-Soap- Oil 


Fig. 1. The interface of two phases (oil, air and water) is an electrical double 
sheath at which charged or polarizable molecules or ions (in this case soaps) 
gather. See W. W. Niven, Fundamentals of detergency, 1950. 
Reinhold Publ. Co., N. York. 


particles increases. If the soap molecules have a length of say 25 A, 
then the particle diameters may increase from 50 A to 100A, giving 
a lowering of the breakdown strength from 25 kV/mm to 9 kV/mm. 
Thus the formation of soaps in the long run will be harmful to the 
insulating oil, even if no resins are present. It all depends on the 
interface energies between the three or two substances (see ref. 3)). 

Van der Minne, Hermanie, Van der Waarden, Mackor 
and Van der Waals jr.4) studied the stabilization of carbon black 
particles by adsorption of soaps or oxidation resins in mineral oil. 
Van der Waarden assumed a mechanism of attraction different 
from the London-Van der Waals forces. He assumed the surface 
area of the carbon black particles to be densely occupied by polar 
carbon-oxygen bonds, which supposition was proved by an analysis 
of the mixture of the adsorbed gases. The product contained 9.2 
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per cent of the unavoidable water and 6 per cent (by weight) of 
oxygen, of which one half could not be removed by heating at 200°C 
and was assumed to be chemically bound. 


S Toya 2S e 
Wee 
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Fig. 2. Different types of micelles formed around an oil droplet in water 
or a water droplet in oil. 


Those dipoles formed by carbon and oxygen attract polar or 
highly polarizable molecules. In the case of a soap molecule the 
inorganic head of the molecule will be preferentially adsorbed to the 
carbon particle (or the water droplet), whereas their organic tails 
may be sticking out tens of A’s (see fig. 2). 
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Aromatic parts of molecules show a far greater interaction with 
polar molecules or particles than aliphatic molecules, which may be 
ascribed to the pronounced polarizability of aromatic nuclei in the 
plane of the rings (Kerr-effect). Thus aromatic parts of molecules 
or metal ends of soap molecules will be strongly adsorbed by carbon 
particles or to water droplets, their aliphatic side chains protruding 
like tentacles from the carbon or water surface. If they are sticking 
out like thorns, the surfaces of carbon particles are being kept apart 
by steric hindrance at a distance of 2/ (J is the length of a side chain) 
and flocculation will hardly be possible if 2 1 exceeds a value of 40 A, 
which is about the diameter of the Van der Waals sphere of 
attraction. 

Adsorption of rods of 20 A length was shown to stabilize a sus- 
pension of carbon particles of about 5000 A diameter, whereas with 
rods of 10 A rapid flocculation occurred (see ref. 4) and 5)). 

The action of aromatic colloid-chemical stabilizers or inhibitors 
in insulating oils may consist in a prevention of flocculation of 
impurities by supplying the impurities with a protective cover 
against further flocculation. It is no wonder that paint and dye 
research has supplied many inhibitors for the construction of cables, 
transformers and capacitors. Both strive to supply a protective 
cover to fibre-like solid hydrocarbons, see § 4. 


§ 3. Electrolyte contents of insulating otl and flocculation of particles. 
Electrolytes may be present in insulating oil and are responsible 
for a certain electrical conductivity different from zero. The work of 
several investigators, as cited by Van der Minne, Van der 
Waarden #2) and Koelmans 5) had made it clear that ions are 
present in liquid hydrocarbons. It is not clear to the author whether 
these ions always consist of ionized oil molecules or whether they 
should be attributed to the presence of foreign matter. 

A number of investigators 4)5) have shown that a soap molecule 
of calcium di-isopropyl salicylate (abbreviated to Ca-‘‘dips’’) in 
oil, if adsorbed to a carbon particle, may be dissociated, imparting 
a positive charge to the particle, whereas tetra-isoamyl ammonium 
picrate (‘‘tiap’’) gives it a negative electric charge. This corresponds 
to the ions adsorbed and the counter ion atmospheres of the particles 
in aqueous suspensions, which may produce attractive or repulsive 
forces. As the carbon-soap and carbon-oxidation resin micelles have 
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opposite charges, flocculation may occur at every collision of these 
micelles. 

Thus, dips may stabilize a carbon-in-mineral oil colloid, tiap may 
do likewise, but if applied together in a certain ratio (say dips to 
tiap in a molar ratio of 20 to 1), immediate flocculation may follow. 
In these experiments tiap was used instead of oxidation resins. 
The latter have a composition which is less well defined than tiap. 

A few (say 1 to 10) millimoles per liter of these additives may 
lead to a stabilization or a flocculation of the suspension in insulating 
oil. If the particles in the oil are kept apart at a distance of 2/ 
(~ 40 A), the weak attractive forces will be influenced by soluble 
soaps and oxidation resins in a very small concentration, which is 
in sharp contrast to the conditions of stabilization and flocculation 
in aqueous suspensions. 

Flocculation leads to a lowering of the breakdownstrength, and 
thus the formation of soaps in the presence of an initial surplus of 
oxidation resins, leading to flocculation, may result in a complete 
breakdown of the insulating quality of the hydrocarbon oil. 

At any rate, addition of a minute trace of dissociated soluble 
matter may easily be the cause of a marked increase in the electrical 
conductivity of the oil. In this respect the measurements of Van 
der Minne and Hermanie 4) and Koelmans 5) on liquid hydro- 
carbons such as mineral oil, benzene and xylene are important. In 
pure xylene the dissociation constant is of the order of magnitude 
10-48, the conductivity < 10-14 ohm! cm-!. Addition of a non- 
ionic stabilizer did not cause the conductivity to increase apprecia- 
bly, but addition of an ionic stabilizer such as ‘‘dips”’ or “‘tiap’”’ of 
a concentration of 10 millimoles per liter increased the conductivity 
to a few times 10-10 ohm~!cm-1, which means an increase by a 
factor 104. 

Van der Minne and Hermanie investigated the stabilization 
of carbon particles in lubricating oil by small amounts of aluminium 
naphtenate, Ca sulphanate, Ca of Mg salts of liquid resin, of phenyl 
stearic acid or substituted salicylic acid. 

Systematic research on fine dispersions in xylene has been 
carried out by Koelmans 5) who used Fe203, AlsOs, SiOz, G 
TiOz, CaCO3, BaSOa, Se and Hg$S particles having a diameter of 
1 micron. Stabilizers applied were non-ionic (oleic acid, stearic acid 
and stearyl alcohol) or ionic (Cu oleate, tri-isoamyl-ammonium 
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picrate (tiap) Ca di-isopropyl salicylate (dips). The concentrations 
were 10 millimoles per liter and in the case of an ionic stabilizer the 
conductivity increased by a factor of 104. 

Van der Minne and Hermanie point to the low ionic con- 
centrations necessary for a 1 : 1 electrolyte to cause or to oppose 
flocculation in a suspension of benzene, which is a factor 10-7 or 
10-8 less than in water. In the latter a concentration of say 20-100 yu 
moles// of a 1 : 1 electrolyte is required for flocculation. In benzene 
400 uw moles// of tiap is required, about 100 times less. Ionic con- 
centrations are much lower because the degree of dissociation of 
the electrolyte is about 10-5 to 10-® as against unity for water. 
Ionic concentrations necessary for a 1:1 electrolyte to cause 
flocculation in benzene are therefore 10-7 to 10-8 times less than 
in water. 

The mechanism will be quite different in the two cases; the 
concept of compression of the double sheath assumed for water, is 
entirely out of the question at the low concentrations occurring in 
benzene. 


§. 4. Inhibitors. If bridges are being formed in the oil, deteriora- 
tion takes place due to local heating in the bridges, by which gas 
bubbles and carbon particles may be formed. In these gas bubbles 
corona discharges take place, by which (in the presence of oxygen) 
acids may be formed. Those corona discharges and the formation 
of carbon may be prevented by keeping the colloid particles small 
and dispersed. Formation of bridges is then impossible. 

It has been shown in §2 that aromatic molecules are surface- 
active towards carbon particles. A number of synthetic aromatic 
inhibitors may be added to supply a protective cover around carbon 
particles, protecting them against further flocculation. It is not 
by chance that dye research has supplied a great number of inhibi- 
tors to the constructors of cables, transformers and capacitors 
(anthraquinone). An extensive review of this subject has been 
compiled by Callinan §). 

Mineral oil may contain a large proportion of aromatic compounds 
possessing long aliphatic side chains. It has become clear that 
purification of mineral oil should not lead to a removal of these 
natural inhibitors, as has been pointed out by Imhof’). The oil 
may deteriorate rapidly if the purification has gone too far. 
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Carbon particles in aromatic liquids are kept dispersed and do 
not flocculate rapidly. The stability of a synthetic aromatic compound 
such as chlorinated diphenyl is clear. 


§5. Liquids free from inhibitors. The author carried out an 
experiment in which any stabilizing influence was assumed to be 
absent. In twice-distilled carbon tetrachloride a carbon sol was 
prepared by effectuating a number of electrical breakdowns. A 
sol of brownish hue ensued. After a few minutes the carbon particles 
in a drop of the liquid on a glass plate were found by microscopic 
observation to grow to larger complexes which even became visible 
to the naked eye. 

In this case no stabilization could take place. Evidently carbon 
tetrachloride should be considered to be less suitable as an insulating 
liquid, although the flow of the coarsening carbon particles through 
the pores of insulating paper may be relatively slow or even com- 
pletely impossible. 


Received 27th October, 1960. 
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MAGNETIC MEASUREMENTS ON PURE MnCO,, 
PURE MnO AND DISSOCIATING MnCoO, 


by J. A. POULIS, C. H. MASSEN and P. VAN DER LEEDEN 


Technological University, Eindhoven, Netherlands 


Summary 


With a sensitive susceptibility balance measurements have been carried 
out on dissociating MnCO3. The magnetic susceptibility of pure MnCOs, 
which was the initial material, proved to follow the Curie-Weiss law; pure 
MnO, as a result of the complete dissociation, showed deviations from this 
law. The recombination of MnO with the collected CO, was studied, the 
results leading to a hypothesis concerning the reaction mechanism of dis- 
sociation and recombination. 


§ 1. Introduction. Measurements are reported of the magnetic 
susceptibilities of pure MnCO3, pure MnO and of the solid products 
of MnCOs, dissociated in different degrees. The ratio of MnCOg 
and MnO in the sample was varied by means of the following 
chemical reactions: 


MeO ewe SANTO EC Ons dee {1} 
MnO + COzg + MnCOs. (II) 


Our aim was to collect data about the magnetic behaviour of 
the pure substances as well as to study the reactions (I) and (II). 

For measuring magnetic susceptibilities according to Faraday’s 
method we used a torsion balance, its deviations being determined 
with an electronic displacement meter. To test the reliability of 
the apparatus during several weeks, a diamagnetic substance was 
mounted in a rather weak magnetic field (3000 Oe). The magnetic 
force measured by our balance was reproduced better than 0.1%. 
Readjustment of the probe caused a deviation of the order of 1%. 
To avoid the influence of such a disturbance of the balance, we 
let the reactions proceed in a container which was permanently 
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attached to the balance during the whole experiment. Therefore the 
original pure MnCOsg was put in a glass bulb, having a connection with 
a much larger one in which the COzg was stored; both bulbs were 
attached to the balance. The first glass bulb could move in a small 
furnace placed between the pole pieces of the magnet while the bulb 
containing the COz was placed outside both furnace and magnetic 
field. When the temperature is raised, reaction I dominates, while 
lowering of the temperature favours reaction II. So we could vary 
the dissociation ratio. 

As neither of the susceptibilities of the pure solids is better 
known than 1% and on the other hand calibration of our balance 
with a better known substance also brought along a comparable 
error, the accuracy of 0.19%, as suggested by the data given, has 
significance only for the relative values of all the data presented 
in this paper. 


B 
ha || 


Fig. 1. A) Vertical en B) horizontal view of the apparatus (schematically). 
For the meaning of the symbols see text. 


§ 2. Apparatus. The apparatus is shown schematically in fig. 1. 
Two crossed bars are suspended from a platinum-nickel strip 4. 
The sample was mounted on the end of the first bar 3, The moving 
part of the displacement meter was attached to the opposite end 
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of this bar. One end of the other bar was connected to a compen- 
sating-force unit, while the remaining end only served for balance. 
The sample was put into an initially evacuated glass bulb 2 moving 
inside a furnace 8, which was mounted between the pole pieces of 
a small electromagnet 9. 

The pole pieces were designed in such a way that the differences 
between values of H @H/éx over the whole area of the field in 
which the sample moved were sufficiently small. The COzg formed 
by the reaction I was gathered in the glass bulb 10 outside the 
furnace and the magnetic field. 

The moving part of the displacement meter consisted of two 
ferroxcube rods 5. These could move in and out of two coils 6. 
The ratio of the coefficients of selfinduction was measured by 
means of an impedance bridge. When this ratio did not vary any 
more, the forces on the balance were taken to be in equilibrium. 
The compensating-force unit consisted of a coil 7 moving in the 
field of a permanent magnet, the current through the coil being 
varied until the induced Lorentz force compensated the magnetic 
force on the sample. 


§ 3. Measurements. Measurements were started on a sample of 
MnCOs3 in the evacuated bulb. The curve of the susceptibility of 
MnCOsg3 versus temperature was measured several times. This curve 
is shown in fig. 2 (curve I). These data and all susceptibility data 
following in this paper have been corrected for the diamagnetism 
of the glass bulb, for the paramagnetism of the surrounding air 
and also for the diamagnetic part of the susceptibility of the sample. 
As long as 370°C was not surpassed, this curve reproduced. As we 
did not expect any equilibria that would adapt swiftly to a dissoci- 
ation degree corresponding to the temperature, we concluded that 
no dissociation had taken place so far. When the temperature 
was raised above 370°C, dissociation started according to reaction I. 
This was concluded from the variations of the measured susceptibi- 
lity at constant temperature. 

When after a dissociation period at a temperature above 370°C 
the sample is kept at temperature below 370°C, recombination 
occurs according to reaction II. This recombination was thoroughly 
studied at voom temperature. In order to obtain information 
concerning the dependence of the susceptibility of mixtures of 
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MnCOz3 and MnO on the temperature during the recombination at 
room temperature, short interruptions were introduced several 
times. During each interruption the temperature was raised up to 
about 300°C and lowered again to room temperature. The sus- 
ceptibility measurements taken while raising the temperature were 
plotted and a curve was drawn through the points. Data gathered 
during the lowering of the temperature corresponded to points 
belonging to the same curve. This led to the conclusion that these 
interruptions did not affect the structure of the mixture, thus 
allowing to determine its y-T-curve. Several subsequent dissoci- 
ations and recombinations were carried out with the same sample. 
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Fig. 2. Inverse susceptibilities versus temperature for MnCOg3 (curve I) and 
MnO (curve II). The lines drawn connect our own measurements. Other 
points are from © Krishnan!), @ Bizette‘), G Birckel?), A Honda). 


The last dissociation period was taken at 500°C. This temperature 
was maintained until the susceptibility remained constant, indi- 
cating that the whole sample was transformed into MnO. This 
complete dissociation was checked by further raising of the tempera- 
ture, no indication of a further dissociation being found in this way. 
Later on we got another check by X-ray analysis. The susceptibility 
of this MnO is plotted as a function of temperature in fig. 2 (curve 
II). These measurements were taken in a similar way as described 
above. Again the curve proved reproducable. This curve too will 
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be discussed later on. With the help of the susceptibility curves 
of pure MnO and pure MnCO3 we are now able to calculate the 
dissociation ratio of our sample at any time of the experiment. 
This can be done at every temperature by presupposing linear 
additivity of the susceptibilities of the components at that tem- 
perature. The complete experiment can now be described with the 
help of fig. 3. In this figure, the ratio of dissociation during the 
whole experiment is plotted against time. 

Several dissociation ratios plotted there were calculated from 
susceptibility measurements at different temperatures between 
room temperature and 370°C without any indication of the tem- 
perature influencing the dissociation ratio during the short time 
of measuring as explained before. The heating periods during 
which dissociation occurred, can be found as discontinuities of 
the curve in fig. 3. 
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Fig. 3. The dissociation ratio versus the time throughout the whole experi- 
ment; points are measured during recombinations at room temperature, 
except points given by the symbol © which were drawn during 
recombination at 250°C. 


The form of the recombination curves indicates that two pro- 
cesses with quite different time constants play a role. The dissoci- 
ation period is always directly followed by a relatively fast re- 
combination of limited extent and duration, and later on by a 
relatively slow recombination. A part of the last recombination 
was carried out at 250°C. 


§ 4, Results and discussion. are 
I. The temperature dependence of the susceptibility 
of MnCOs3. According to the Curie-Weiss-law the 1/y — T curve 
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for a paramagnetic substance should be a straight line, given by 


14.052 
= ee 0), 
16 Py? 


where Py is the magnetic moment of the paramagnetic ion expressed 
in Weiss magnetons and 6 a constant characteristic for the sub- 
stance. 

In fig. 2 (curve I) the results of our measurements after cor- 
rections mentioned previously in this paper are plotted. Our graph 
shows a straight line in the case of MnCO3. From the separate 
points we calculated the constants P,», using the method of least 
squares: 

Py = 28.74 £0.09, @ = — 71.4 + 3.9°K. 


The given uncertainties are estimated standard errors. For com- 
parison we give the results of resp. Krishnan 1) and of Birckel ?) 
for MnCOs3: 


Krishnan: * 2a = 2/.0'-- Le. 6 = — 13° + 60°K; 
Birckel: Pe 250 0:25. 6 = — 56° + 9°K. 


Significant incompatibilities do not occur. 


II. The temperature dependence of the susceptibility 
of MnO. In fig. 2 (curve II) our measurements concerning MnO 
are given and compared with those of Honda 8) and of Birckel 2). 
The deviations from the Curie-Weiss-law, suggested by the results 
of Honda, became doubtful after the work of Birckel; they 
are however significant according to our measurements. Recently 
results concerning the susceptibility of MnS have been published 5). 
Here similar deviations from the Curie—Weiss-law are stated. 

We tried to describe our results by putting 1/y equal to a 
quadratic polynomial in T and determining the three constants 
by the method of least squares. The discrepancy between the curve 
thus found and the measured points could be described by a stan- 
dard error of 0.08%. Deviation from a straight line would be 
systematic and amount to about 5%. A study of these discrepancies 
taught that higher than second order approximations could not be 
expected to give a significantly better description. 

Values of Py were calculated at various temperatures by using 
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the Curie-Weiss-law for the tangent (fig. 4). The values of P,, 
plotted in fig. 4 below room temperature were the result of extra- 
polation of our polynomial. They show the quite remarkable fact 
that near the Néel temperature the number of magnetons of MnO 
is equal to the theoretical number P,, — 29.3 corresponding to 
ies /2. 
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Fig. 4. The temperature dependence of the number of Weiss magnetons 
of MnO versus temperature. The line drawn connects our own measurements, 
the dotted line is an extrapolation. The other points are measured by 
a Birckel®)) A Bizette,4): 


Ill. The reactions. From our magnetic measurements in 
conjunction with X-ray and microscopic investigations some 
conclusions about the reaction mechanisms may be drawn. 

Microscopical observations showed that the grain sizes were the 
same before and after the measurements described in fig. 3. The 
X-ray investigations showed that we had to do with large and small 
crystals of both MnO and MnCOs. The large crystals (> 500 A°) 
gave sharp diffraction lines while the small crystals (~ 300 A®) 
gave broader ones, which were displaced a bit, due to the lattice 
constants being influenced by the size of the crystals. X-ray 
investigations have been made before and after the measurements 
of fig. 3. At several intermediate points X-ray information was 
drawn from a similar sample, which was given the same heat 


treatment. 
The X-ray experiments showed that the original MnCOg con- 
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sisted of large crystals. After each heating period 9% of the total 
amount of Mn ions formed small MnO crystals. The fast part of a 
recombination resulted in the formation of 9% of small crystals 
of MnCOg3 and elimination of the small MnO crystals. 

This led us to the following suggestion about the reaction 
mechanism: An original MnCOg grain consists of large crystals. 
By heating, the surface layer the grain can easily be dissociated ; 
small MnO crystals are formed. The central part of the grain can 
only be affected by more extensive heating. 

The small MnO crystals in the surface layer when in contact 
with CO, gas at room temperature can easily recombine to form 
small crystals of MnCOxg, the MnO crystals at the inside of the 
grain recombine much more slowly. 
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COUETTE FLOW OF A FULLY IONIZED GAS, 
CONSIDERED AS A TWO-COMPONENT FLUID 


by L. A. PELETIER *) and L. VAN WIJNGAARDEN 


Mededeling No. 94 uit het Laboratorium voor Aero-en Hydrodynamica van de 
Technische Hogeschool te Delft, Netherlands 


Summary 


The equations, governing the behaviour of a fully ionized gas, as given 
by Spitzer, are applied to two types of magnetohydrodynamic Couette 
flow. The features of the flow are expressed in terms of the Hartmann 
number and a parameter g, being the ratio between Larmor frequency 
and collision frequency. Compared with the results of the one-component 
theory an additional velocity component is found. 


§ 1. Introduction. In dealing with problems concerning the motion 
of an ionized gas through a magnetic field one may consider the 
gas to be continuous as long as the mean collision frequency is large 
with respect to the gyration frequency of the motion of charged 
particles about the lines of magnetic induction. 

When the ratio between gyration frequency and collision fre- 
quency is not negligibly small, additional effects must be taken 
into account. This is done, following Spitzer 1) by considering a 
fully ionized gas as a mixture of two components, one consisting 
of positively and the other of negatively charged particles. From 
the equations of motion of the two components the equations 
governing the macroscopic velocity and current density of the gas 
can be deduced. 

In § 2 we apply Spitzer’s equations, extended with a viscous term, 
to the flow of a viscous fully ionized gas between two parallel 
infinite plates, one being at rest and the other moving with a 
constant velocity. A constant magnetic field is applied in the 
direction normal to the plates. We consider two different cases of 


*) At present at the Massachusetts Institute of Technology, Cambridge Mass. U.S.A. 
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this Couette-flow. In § 3 the fluid is shortcircuited parallel to the 
plates whilst in § 4 both plates are insulators. 

A velocity component, not present in the one-component theory, 
perpendicular to the normal to the plates and the velocity of the 
moving plate is found. An additional shear force in that direction 
results from this velocity component. This effect is shown to be 
related to that noticed by Broer and the present authors 2) in 
studying the forced motion of a non-viscous fully ionized gas between 
two parallel plates at rest. 


§ 2. Equations for Couette-flow in two-component theory. Consider 
two parallel plates of infinite length and width. One of them coincides 
with the x*z* plane, is at rest and is a perfect conductor. The other 
one, a non-magnetic insulator, moves at y* = A with a constant 
velocity U in the direction of the positive x* axis. Outside the plates 
is a constant homogeneous magnetic field directed along the 
positive y* axis. 

The space between the plates is occupied by a fully ionized 
viscous gas, consisting of protons and electrons. Moreover we assume 
the gas as a whole to be neutral and in thermal equilibrium. It is 
assumed that U is low enough to make compressibility effects 
negligible. We seek the velocity distribution in the gas and the 
forces on the plates, when the equations of Spitzer’s two-component 
theory are applied to the motion engendered by the moving plate. 

First we have to consider that the equations given by Spitzer 
(ref.1), p. 21) result from the linearized equations of motion of 
the components, and that they give no account of viscosity. The 
first point gives no difficulty because the convective acceleration, 
neglected by Spitzer, vanishes in our case as all quantities depend 
only on y*. 

As for viscosity, Spitzer’s equations can be extended to viscous 
problems by introducing in the equation of motion the viscous 
term of the Navier-Stokes equation. Doing this we have, bearing 
in mind that the motion is steady (in M.K.S. units): 


Jx B—VP+7V%=0 (1) 
and 


a | J <°B VP 
E*+ yx B— — —- —., (2) 
oO 
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where 

= velocity 

= magnetic induction 

pressure 

= electric field strength 

= charge of proton 

= number of particles of each component in a m3 
electrical conductivity 

current density 

coefficient of viscosity. 


Subs 
Il 


* 
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(1) and (2) differ from the corresponding equations in the one- 
component theory by the terms on the right-hand side of (2). They 
represent the influence of the Hall-current and the partial electron 
pressure respectively. In the one component theory the right hand 
side of (2) is zero giving Ohm’s law. It can be made plausible that 
in (2) viscous terms can be neglected, the electronic contribution to 
viscosity being very small (cf. Chapman 3)). 
In addition to (1) and (2) we have 


Y-ye= 0; (3) 
Vv -B=,0, (4) 
V3 El), (5) 
VY Xx By, (6) 


where yu is the magnetic permeability. In order to formulate the 
problem in terms of dimensionless quantities we take 


x* — xh, y* = yh, z* = zh,v = uu, B = bBo, 


; nU 
E* = EUBp, J = joUBo, P=p——. 


Then (1)—(6) become 


MY x b—Vp+V2u =0, (7) 
Beeb) 0 =f = 0, (8) 
Boa =40. (9) 
bea, (10) 
ix E= 0, (11) 


Vx b= Rnj (12) 
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These equations contain three dimensionless groups: 
$ 


M = Boh ) Se eun 
7" 
oBo 


en 


(13) 


Of = 


M is the Hartmann number, R» the magnetic Reynolds number, 
while g is the ratio mentioned in the introduction. This is easily 
seen by remembering that kinetic theory gives o = e2n7/m, where 
7 is the mean collision time for collisions between electrons and 
protons and m is the mass of an electron. Upon inserting this in 
(13) we obtain ¢g= 7, where w = eBo/m is the frequency of 
gyration about the lines of magnetic force for an electron. 

The quantity g measures the number of rotations about lines of 
magnetic induction which an electron performs between collisions 
with protons. Here we must restrict ourselves to small values of ¢ 
as (1) and (2) are valid only for this case. For larger values of g the 
quantities o and 7 cannot be represented by scalars. 

Next we have to formulate the boundary conditions. They must 
state that there is no slip at the plates, that no electric current is 
flowing into the upper plate, this being an insulator, that the 
y-component of the magnetic induction is continuous at both 
plates and that the other components are continuous at the upper 
plate since this one carries no current. 

Hence we have 

y == 0: = 0,-by = I, 
Yl: tg = 1) Uy SS Uz =0, fy =O, Oy = 15 bp Hb, = 0. 
From the boundary conditions and (9)—(12) we deduce, keeping in 
mind that there is only dependance on y: 


Uy = 0, by = 1, J7y = 0, Ez = constant, FE, = constant. 
Using these results, the x and z components of (7) become 


: 224, 
— M4, 4+ Te = 0, (14) 
Mic thee 
Ja dy? , (15) 
and of (8) 
Ex — Uz + Gz — Ja = 0, (16) 


Ez+ Ug — iE = Va = 0. (17) 
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Solution of (14)-(17) requires the knowledge of six constants, 
including Ez and Ez. The boundary conditions for u give us four 
relations between these constants. Hence two additional relations 
are needed. In the next two paragraphs two different choices for 
these relations are made. 


§ 3. Shortcircuited gas. We suppose the gas to be shortcircuited 
in the x and z direction. This implies 


Beet Oh (18) 


It is possible to solve (14)-(17) together with the boundary condi- 
tions and (18) for all values of g. But; as already indicated, the 
validity of the present theory is restricted to small values of gq. 
Therefore we give here the solution for small g. When neglecting 
terms of second order in g, we obtain from (16) and (17) with help 
of (18) 


Uz = Quz — Ix (19) 
or with (15) 
ie 1 du, 20 
Uz = gu Ss - : 
Zz Guz We dy2 ( ) 


The first term at the right-hand side of (20) is related to a result 
obtained by Broer and the present authors 2). They studied the 
motion of a non-viscous fully ionized gas between two infinite 
parallel plates, a constant magnetic field being applied normal to 
the plates. In their problem the motion is caused by a constant 
volume-force. It was shown that the gas moves at an angle y with 
respect to the direction of the driving force where @ is given 
by tg » = 3g. In the present problem the driving force is the x- 
component of the viscous stress. The first term at the right-hand 
side of (20) shows that a similar effect is present here. The second 
term represents the influence of the shear stress in the z-direction 
on the drift velocity uz. From (19) we deduce that 7, is of the first 
order in g. Hence we obtain from (17) and (18) that in our approxi- 
mation uz = 7. Inserting this in (14), yields 


2 
dz — M2u, = 0. (21) 
dy? 


: Obviously there is no influence of g on wz up to the second order. 
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The solution of (21) satisfying the conditions 


y= 024, = 0, y= 1li-ug= 1 


sinh My 
i Lotta oe 22 
a ENE (22) 


With the help of (22) and the boundary conditions 
Ayes Occ a= ol) Siesta 5 ethan) 
the solution of (20) is 
= ARE (cosh M sinh My — y cosh My sinh M). (23) 
2 sinh?M 


For large M, ie. 
cosh Mw sinh MS MS 1, (24) 


we Can approximate (22) and (23) for y > 0 by 
Uy, © exp[— M(1 — y)], (25) 
uz ~ 3Mq(1 — y) exp[— M(1 — y)]. (26) 


These expressions show that for large M the motion is restricted 
to a boundary layer on the upper plate (whose thickness is of order 
1/M), the ratio wz/a in this boundary layer being given by 


Uz/U, ~ 4Mq(1 — y). 


Once wz and wz are known, j can be calculated from (14) and (15), 
b from (12), Ey and df/dy from the y-components of (7) and (8). 
As these quantities are of no great interest, we don’t give them here. 
We only mention that the two integration constants involved in 
the calculation of by and b, from (12) are determined by the con- 
ditions that 6; and bz must be zero at the upper plate. The values 
of 6, and 6, at the lower plate then determine the current in the 
lower (perfect conducting) plate. This current is equal and opposite 
to the total current in the gas. 

Quantities of interest are the viscous stresses exerted at the 
plates. They are given in dimensionless form by the first derivatives 
at the plates of wz and wp. Denoting these dimensionless stresses 
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by kz and kz we have, from (22) and (23) at the lower plate, 


(SS) 
Re —= SSS 
dy: /y=-9 sinh M ~ 


ag ( = Bete (MM pena si 

: dy /y-9  2sinh2M (M cosh M — sinh M), 

and at the upper plate, 
diz 

4 =—( ) = — McothM, 

x dy y=1 (27) 
Hee ( a ) Jot inh 2M — 2M 

Ss dy Jy-1 4 sinh2M isi ri ). (28) 


In the case where M is large (24), (27) and (28) yield for the ratio 
between kz and kz at the upper plate 

Re[Ra = — 39. (29) 
The expression (29) resembles the expression obtained by Broer 
et al., in their already cited note, for the mechanical Hall-effect, i.e. 
the pressure gradient build up perpendicular to the driving force, 
when in their problem the gas is compelled, by walls inserted in 
the flow, to move parallel to the driving force. This pressure 
gradient was found to be 4g times the (volume) force. 


§ 4. Insulated plates. We return to equations (14)—(17) and assume 


now that both plates are insulators, electrically and magnetically. 
Hence the net current in the gas is zero: 


1 1 
Siz dy = 0, fjzdy = 0. 
0 0 


With these conditions we obtain from (14) and (15) upon integration 


between O and 1: 
eye ee 
dy /y=0 dy /y=1 


are 
dy /y=0 \ dy /y=1 

The solution of (14)-(17) then has to satisfy (30), (31) and the 
conditions for wu given in § 3. 
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Here again we consider only small values of g. Neglecting terms 
of second order in g, we obtain, using the same arguments as in the 
foregoing paragraph, 


1 du, 
Uz = Qua + V2 dy? oe Ex Sle gE z (32) 
and 
1 2 
a ce ae ioe ES (33) 
M2 dy? 


The solution of (33), satisfying 


and (30) is 
sinh M(y — 4) ] 
=1/1 : 34 
eee [ Feo dain bag (34). 
E SS = 4. (35) 


Substituting (34) and (35) into (32), we find using (31) and the 
conditions 


Vie Othe 0 nigel goes 
the following expression for wz: 


Mq 


Le ae aL FETE OEE 
4 sinh? 4M 


‘ [cosh $M sinh M(y— $)—(y—}) sinh 4M cosh M(y — 4)], (36) 


whilst Ey = 0. 
When we move the x*-axis to y* = 44 in the dimensional coordinate 
system, then the ordinate is y*’ = y* — th. Upon choosing 


we have 
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Inserting this into (34) and (36) yields 


sinh 4MC 
ay emt eee (37) 
sinh $M 


Ug — 


nie 


Mq 


eg = 
8 sinh? 4M 


[cosh $M sinh 4M¢ — ¢ sinh $M cosh $M]. (38) 
Comparison of these expressions with (22) and (23) shows that in 
the region 0 < € < 1 the expressions for (wz — 4) and uw; as given 
by (37) and (38) are the same as the corresponding ones obtained 
in the foregoing paragraph, provided we assume in § 3 the distance 
between the plates to be $/ instead of / and the velocity of the upper 
plate 4U instead of U. 

This similarity and the antisymmetry of the right-hand sides 
of (37) and (38) with respect to ¢ = 0, enables us to obtain all 
quantities from the corresponding ones in § 3 by applying the scaling 
rules mentioned above. To obtain the stresses on the plates we 
have only to replace M by $M in (27) and (28). Then we have from 
(27) and (28) in view of the antisymmetry 


kz = + 4M coth 3M, 


Mq ( M — ae 
8 sinh2 1M /’ 


keg = 


where + refers to the lower and — to the upper plate. Clearly, 
there is no net force, in agreement with the fact, that the net current 
vanishes. 

When M is large there is a boundary layer on each plate; (37) 
and (38) together with (24) give for the velocity distribution in the 
boundary layers 


Uz —4 ~hkexp[— 34M (42+ 1), 
Mq 


3 (+ 6+ 1) exp[— 3M(4 ¢ + 1)], 


Uz © 
where again + refers to the lower and — to the upper plate. In 
§ 3 we obtained the expression (29) for the ratio between kz and 
kz at the upper plate when M is large. (29) does not contain M and 
therefore holds here too. Hence we have at both plates kz/kz = — 4q 
when M is large. 
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ELECTROSTATIC DIPOLE MOMENT OF A 
DIELECTRIC. CUBE 


by T. W. EDWARDS and J. VAN BLADEL 


Department of Electrical Engineering, 
University of Wisconsin, Madison, Wisconsin, U.S.A. 


Summary 


The electrostatic dipole moment of a dielectric cube introduced into a 
uniform electrostatic field can be determined by solving an appropriate 
surface integral equation. High-speed digital computers have been employed 
to solve the matrix approximation to this equation. 


§ 1. Introduction. The dipole moment of a dielectric body 
introduced into a known electrostatic field can be calculated by 
separation of variables for a few simple shapes such as the cylinder, 
the sphere, and the ellipsoid. However, there exists a surface integral 
equation method of calculating the dipole moment which is valid 
for all shapes. The actual solution of this integral equation must 
be effected by approximate methods of considerable complexity. 
The advent of modern high-speed digital computers permits solution 
of many problems previously intractable. The authors decided to 
apply these methods to the solution of the dipole problem. For the 
sake of simplicity, a cubical shape in a uniform field was chosen. 
The results obtained are applicable to the scattered field of a 
dielectric cube immersed in a plane wave whose wavelength is large 
compared with the dimensions of the cube. 


§2. The integral equation. Phillips!) has shown that the 
electrostatic potential distribution ¢ inside and outside a dielectric 
body inserted into a pre-existing potential distribution y (see fig. 1) 
can be expressed, respectively, in terms of the resulting surface 
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potential distribution by the formulae 


—K 0 1 
es(P) =P) + | | 6 kt 


An any |Fp — Fp, 


—K 0 1 
RS aS ee 
An [Jo Ong rp ep eee 2) 


where « is the dielectric constant of the body. 


1 
$(Pe) = p(Pe) + 


Fig. 1. A dielectric body bounded by the surface S. 


Allowing the points P; and Pe to approach a point Pg on the 
surface S of the dielectric body leads to the following integral 
equation for the surface distribution: 


k+ 1 k— | 0 1 
vib) = “2 40) +=] 4) = as, @) 


§ 3. Specialization to a dielectric cube immersed in a uniform 
electrostatic field. The dipole moment of the cube obviously 
depends on the orientation of the cube with respect to the direction 
of the initial uniform field Eo. However, by separating the field 
into its three components, each normal to a pair of faces of the 
cube, it is seen that the problem reduces to the calculation of the 
dipole moment induced by any one of the three components. 

Thus, considering the z orientation only (see fig. 2), the original 
potential distribution yp is given by 


y= Eoz. (4) 

The surface integrals in (1), (2) and (3) may be approximated 
by finite sums by dividing each face of the cube, of total area 4a2, 
into N? equal sections, and by considering the potentials 4 and y 
to be constant over each section. In this manner (3) gives rise to 
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a matrix equation 
k+ 1 l—k 
(v(bo)] = “> — [#(G0)] + — [Alis(o, (6) 


where [A] is a square matrix of order Q. Consideration of the 
symmetry properties of the surface potential distribution 4(f) 
reduces the order of Q and the number of unique values of potential 


Each Side Divided Regions of Unique Potentials 
Into N* Elements z  t(XY+) (3) 
Elements Elements 


Fig. 2. A dielectric cube immersed in an otherwise uniform 
electrostatic field. 


from 6N? to (when N is even) 3N2/8 + N/4. The storage capacity of 
the computer employed (an IBM 650) limited the maximum usable 
value of N to 10, i.e. Q = 40*). Obviously, the matrix approximation 
to the integral equation improves in accuracy as WN increases. 


Fig. 3. Total potential distribution (along the z axis) of a dielectric cube. 
Note that, in this and succeeding figures, the properties of the cube approach 
those of a metallic cube, as the dielectric constant « increases indefinitely. 


*) For more details on the establishment of (5) and for details of the computational 


program, see 2). 
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Fig.3 shows the total potential distribution along an axis 
through the centre of the cube and extending in the direction of 
the initial field. Fig. 4 shows the fraction of this potential distri- 
bution due to the presence of the dielectric cube. 


10 


oA 
$ 
aE, 
r x 
0.01 h 
loo 
20 
10 
t 3 
4 
2 
Zz Li 42 WS 
Ee err rae ae 
1 2 Bonk 7 2105 7 *20 


Fig. 4. Dipole potential distribution (along the z axis) of a dielectric cube. 


§ 4. Dipole moment. The dipole moment of the dielectric cube 
can be read directly from fig. 4, where the 1/R2 decrease of the 
potential at large distance is evident. It can also be determined 
from the definition of the dipole moment P and from (2) as 


peer Be, ah P 6 
PDipole = Aare es re [p(Pe) — (Pe)].- (6) 

This yields 
P= en(« — 1) ffs $(b) Ap ASp. (7) 


For our configuration the dipole moment will be in a direction 
opposite to that of the initial field. Its magnitude is shown in 
fig. 5, where a comparison of the computed dipole moment of a 
dielectric cube and those of two significant spheres, the inscribed 
sphere and the equivolume sphere, is evidenced. 
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Fig. 6 shows a comparison of the dipole moments per unit volume 
of a sphere and of a cube, and gives an indication of the effectiveness 
of these two shapes in producing dipole moment. 
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Fig. 5. Dipole moments of dielectric bodies. 
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Fig. 6. Dipole" moments (per unit volume) of dielectric bodies. 
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SKEW RAY TRACING THROUGH TORIC 
REFRACTING SURFACES 


by J. BECKER 


N.V. Optische Industrie De Oude Delft, Delft, the Netherlands 


Summary 


A simplified method for skew ray tracing through toric refracting surfaces 
is described and formulae are given. 


§ 1. Introduction. In an optical system having only concentric 
refracting and reflecting spherical surfaces all aberrations due to 
lack of symmetry are absent. A system of only three components 
may constitute a wide angle objective of reasonably large aper- 
ture !). The image plane is spherical too, and for an infinite object 
it is concentric with the surfaces of the optical system. Used as a 
camera objective, the film should adapt to this spherical focal 
plane. With sufficient accuracy this is only possible for a relatively 
small field of view. However, the image surface can be made 
cylindrical by introducing a field-flattening lens. In that case the 
accurate positioning of the film presents no difficulties. When 
flattening the image plane for a vertical section (fig. 1), the field 
of view in the horizontal plane still has no limitations, whereas 
now the useful vertical field of view is limited by the unsharpness, 
which is due to the aberrations introduced by refraction through 
the field-flattening lens. This lens has a toric surface and a cylindrical 
surface on which we want the image to be formed. A method for 
ray-tracing through toric refracting surfaces has been described 
by Murray 2). As will be shown, a simplified calculation scheme 
may also be used. 


§ 2. Determination of the point of intersection. We first determine 
the point of intersection Q of a ray, coming from a point P on the 


— 156 — 


SKEW RAY TRACING THROUGH TORIC REFRACTING SURFACES 157 


preceding surface of the system, with the toric surface. At the 
top T (fig. 1) the toric refracting surface has radii of curvature 7 
and 7;¢. The distance between the centres of curvature a = 7, — 7. 


cylindrical 
imoge plane 


Fig. 1. Cross-section of the toric refracting system. 


The system of coordinates has been taken such that the origin O 
coincides with the centre of curvature with the radius 7. The x-axis 
coincides with the optical axis, the y-axis is the axis of the cylindrical 
image surface. The toric surface is given by the equation 


f(x, y, 2) = (x2 + y2 + 22 + a? — 7?)2 — 4a2(x? + 22) = 0. (1) 


For the incident ray through the point P (&, 7, ¢), with cosines 
L, M, N relative to the x-, y- and z-axis respectively, the equation 
reads 


%=&—Lp, y= —Mp, z=£—Np, (2) 


where p is the distance between the point P and the point Q (x, y, 2). 
We find after substituting (2) in (1) and some rearrangements 
that 


p! — Abp? + 2p2(2b + c2 — 2a%d) + 4p(2a%e — bc) + o? — 4a2g = 0, (3) 
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in which 
b=éL + 7M + ON, 


d = L2 + N2, (4) 
é=éL O(N, 
g= e+ C2 


Now p is the skew distance between the two refractive surfaces. 
For the principal ray of the pencil to which PQ belongs, the distance 
may be calculated from a conventional meridional ray tracing. 
Supposing this distance to be ¢, we find, putting p =? + A 


(¢ + A)* — 4b(¢ + A)3 + 
2(¢ + A)2(2b2 + c — 2a2d) + 4(¢ + A) (2a%e — bc) + c2 — 4a2g = 0. (5) 


If A is small, 44 and A? may be neglected, while 4 can be found 
from the resulting quadratic equation, giving a first approximation 
for p. The value of p thus found is corrected in one or two steps by 
insertion into (3), until p has the accuracy required. The coordinates 
of the point of intersection QO of the ray and the toric surface are 
then given by (2). 


§ 3. The refraction. The cosines of the angles with the co- 
ordinate axes for the normal in Q (x,y,z) on the toric surface 
are given by 

Met 1 & ise 


woe’ way’ wa’ ©) 


Palen ianles) 
ee ay) Tae) 
The angle of incidence » between the ray and the normal in Q 
is found from 


MEN Kumeecar ue. 


in which 


The angle of refraction y’ is obtained by application of Snell’s 
law. Hence, the value of C = n’ cos ’ — n cos p is determinable. 
This value is needed for the calculation of the direction of the 
refracted ray. 
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The cosines L’, M’ and N’ for the refracted ray in the coordinate 
system of fig. 1 may be determined as follows (fig. 2). Take a length 
of line QA = m on the incident ray, a length of line OB = n’ on 
the refracted ray. Then AB is parallel to the normal (fig. 2). If 


Fig. 2. Determination of the cosines for the ray, 
refracted at the toric surface. 


the broken line 0A~AB-BQ is projected on to the coordinate axes, 
the sum of the projections on the x-, y- and z-axis is in each 
case equal to zero. The cosines and lengths of the three parts 
of the broken line are given in table I. 


TABLE I 
| cosines | lengths 
QA L M N n 
ier weelrhore Iie Roy 
AB Paes = a aon n’ cos yp’ — ncosp = C 
BQ iZ M’ N’ | x’ 


howe 
Pile ies Sep Tio 
OX 
or 
(Gam lac 
pp sah Shi @) 
n n' w OX 


Projection on to y- and z-axis respectively gives 


Coot da 6 

M’ = et M + 7 / ? (9) 
n n’ w oy 
n Cook cf 


SNS (10) 
nN 
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$4. Conclusion. In this way the refracted ray is completely 
determined, and its intersection with the image plane can be 
found. In practice with a desk computer the ray tracing of seven 
rays from one pencil was achieved in three hours, making it tairlv 
easy to obtain the best correction. The formulae are also suited 
for programming the calculating scheme for an electronic digital 
computer. As a direct application the design of a new(horizon)- 
camera *) having a horizontal field of 360° and a cylindrical image 
plane may be mentioned. 
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ON THEOREMS OF MINIMUM ENERGY 
DISSIPATION IN MAGNETOHYDRODYNAMICS 


by L. N. TAO 


Department of Mechanics, Illinois Institute of Technology, Chicago, Illinois, U.S.A. 


Summary 


This paper is concerned with some theorems on the dissipation of energy 
in magnetohydrodynamics. It is shown that when certain conditions, (2.21) 
or (2.27) of the text, are satisfied, the steady motion of an electrically 
conducting incompressible fluid has an absolute minimum of energy dissi- 
pation. Furthermore, when the same conditions are satisfied, an unsteady 
motion with steady boundary conditions always tends to its steady state, 
which is stable as well as unique. The present theorems are also applicable 
to ordinary hydrodynamics and magnetohydrostatics. A brief implication 
of these theorems is discussed. 


§ 1. Introduction. In studies of magnetohydrodynamics many 
authors have considered the situation of equipartition of kinetic 
and magnetic energy. The suggestion is based on the physical 
argument that fluid motion tends to stretch the magnetic lines of 
force, and so to increase the magnetic energy. The magnetic 
energy, however, cannot be increased indefinitely. When the 
magnetic energy becomes too large, the tension along the lines of 
force checks any further extension of them. Thus in an ultimate 
steady state the inertial forces can be expected to be comparable 
with the effects of tension along the lines of forces. In other words, 
an equipartition between kinetic and magnetic energies or its 
approximate form will finally be reached in a steady state. This 
was first suggested by Walén1). Elsasser?) also commented 
that the equipartition may be conjectured from the symmetric 
form of the governing equations of magnetohydrodynamics. In 
the following we intend to give some theorems which may be 
helpful in discussions on equipartition of energy. 
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In ordinary hydrodynamics, i.e., in the absence of magnetic fields, 
some general theorems on minimum energy dissipation have been 
discussed by Helmholtz 3), Korteweg 4) and Rayleigh®). They 
are generalized in this communication to magnetohydrodynamics. 
In the first part it will be shown that under certain conditions, 
which include equipartition of energy as a special case, the total 
energy dissipation of a motion within a closed system is a minimum 
when compared to any other motion subject to identical boundary 
conditions. Hence a motion with equipartition of energy possesses 
an absolute minimum of energy dissipation. Furthermore, we will 
also show that when the velocity and magnetic fields are maintained 
steady over the boundary, an unsteady motion satisfying previous 
conditions always tends to its steady state. The type of steady 
motion finally attained is stable, as well as unique. The motion 
with equipartition of energy is one of these motions, and therefore 
it is a stable equilibrium. 


§ 2. Steady motion. Consider a steady motion M of an incom- 
pressible fluid in a region D bounded by a closed surface &. Let 
its velocities be w; and magnetic field intensities H;. Then the 
total dissipation of energy of the motion *) is 


Q mes (Gm + $e) dv = f (Tijui,g + 7?/0) do, (2.1) 


D 


where ¢m and ¢¢ are, respectively, the mechanical (viscous) and 
electromagnetic dissipation per unit volume. Tj is the viscous 
stress tensor, 7; the current density **) and o the electric con- 
ductivity. 

Now let us consider another motion M* characterized by velo- 
city u;* and magnetic field intensity H;* which satisfy the same 
boundary conditions as the motion M. Writing u;* = uj + u;' 
and H;* = H; + Hj’, we then have wu,’ and H,’ vanishing at all 
points on the boundary. The dissipation of this alternative motion 
M* is 


*) The summation convention is used in this paper. 
**) Within the validity of magnetohydrodynamic approximation, i.e. negligible dis- 


placement current, the current density may be interpreted with reference to either a 
moving or stationary observer. 
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es { ea sra es 


D 
| Le rere 

= |e + T"%y) (Wig + u'4,3) dv + = (94 + 9a’) (7e + 7:/) dv = 
D D ; 

¥ | oo sii | ges Ca sire 
D D 

i | 

+ | (e+ de) dv + 2] jai ae, (2.2) 

D D 


From the viscous stress tensor 
Tig = (tig + 3,4); (2.3) 
where 7 is the viscosity of the fluid, we find 
Tigw'a,g = n(uig + Uy a)'4g = 
= 4(Uw's,g + U'5,4)Ma,g = Tyas. (2.4) 
The substitution of (2.4) into (2.2) gives 
O*=Q+Q'+ 2 Ty's dv + Ale 1474 dv. (2.5) 
D D 


Now we intend to show that these two integrals vanish under 
certain conditions. Using integration by parts, the first integral 
of the right hand side of (2.5) is 


J Tiwi, W = Titi yay — [mT y,3 do = — f nei'ua, jy dv, (2.6) 
D D DD 
since the u;’ vanish on the boundary and uj; = 0. 


From the momentum equation of an electrically conducting 
liquid in the steady state 


puguig = pPi + Li — fy + Mi, 9, (2.7) 


where p is the density, F; the extraneous force, Lj = weijnyjHe 
the Lorentz force and w the magnetic permeability, we may write 


ST iju'4,5 dv =f ud [pugui,g — HH Hi,j — pli + 0,1 + w(AyA)) | dv. (2.8) 
D 13) 
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Now, we observe that if pF; has a potential 2, 
ply = — Qi (2.9) 


and 


UjUi,g — ms A;Hy,; = yilp, (2.10) 


where y is a single-valued function in the domain, then the integral 
ST yt 63 du 0 (231 1) 
D 
by virtue of 
Jug (o + y) dv = ui'(@ + Y)loay —f(o+y)uv44dv=0, (2.12) 
D D 
o=p)+ wHjAH; + Q, (Zak) 


due to “;’ = 0 on the boundary and w’;,; = 0. 
Using the vorticity vector 6; and Maxwell’s equation between 
7, and H; and taking the curl of (2.10), we then have the condition 


Vx#E x u) = 72 (43 x i) (2.14) 


in place of (2.10). 
The second integral in (2.5) may be written as 


Tag Wes Pali 
[Pie dv =|-H ag(Hi,3 — H; i) dv= — ie Ay Ay dv. (2.15) 
D D D 


Here we have used integration by parts and Maxwell’s equation 
H;,; = 0. The elimination of the electric field intensities E; from 
Maxwell’s equations (steady case) and Ohm’s law 


eqgrEn,3 = 0, 
eige ej = Ji, (2.16) 
Ji = O(Ei + peageujHy) 
gives 2) 
Uj;H;y3 — BLE a — Hy = 0. (2.17) 
ou 
Therefore, if 
ujHi,j — Hywi,g = 7,4, (2.18) 
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where y is another single-valued function in the domain, then 
(2.15) is 


l ae / 
ie Jiji ca Eee [um (ujH iy, T= Hju;,3) dv = [atten dve=s0: (2.19) 
D D D 


We may also replace (2.18) by taking the curl of (2.18); this gives 


Vat x Hy = 0. (2.20) 
Hence when 
pF = — VQ, 
rg x uy 72 (“jx H), (2.21) 
p 
V4(u x A) = 0, 
(2.5) is simply 
er =040 (2.22) 


Furthermore, Q’ is non-negative. It vanishes only when w’;,; = 0 
and H’;;=0. This implies that w;’ and H,’ remain constant 
throughout the domain as if the system moved as a rigid body 
in a “‘frozen’’ magnetic field. Since we have wu,’ and H;’ vanishing 
on the boundary, no such motion is possible. Therefore, the motion 
M which satisfies conditions (2.21) is unique and possesses an 
absolute minimum of dissipation of energy. In conclusion, we have 

Theorem IA. In magnetohydrodynamics (without displacement 
current) the motion of an electrically conducting viscous incompressible 
fluid has a smaller dissipation of energy than any other motion which 
obeys the same boundary conditions, provided conditions (2.21) are 
satisfied. 

In the case of equipartition of energy §), i.e. 


= 
u=/4u, (2.23) 


the motion clearly satisfies (2.21), and hence its energy dissipation 
is an absolute minimum. | 

In theorem IA we have used the fact that the extraneous force 
pF; has a potential Q. This condition may be removed. In that 
case from (2.6) we have 


[ Tiyw's,g dv = — f gus'i,yy Av. (2.24) 
D D 
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The integral vanishes when . 
U4j5 = O (2.25) 
or 
V26 == 0; (2.26) 
This condition is satisfied by Stokes’ flow or any other “‘slow”’ 
motion, i.e., motion void of inertial force effect. 
Writing the equivalence of V2(u x H) =O in the form of 
V2j = 0, we have 
Theorem IB. In the case of the non-existence of a potential Q 
of the extraneous force, theorem IA 1s true, provided the motion 
satisfies 


Ve=0, R=. (2.27) 


§ 3. Unsteady motion. Consider an unsteady motion of an 
electrically conducting liquid in a region D bounded by a surface &. 
The boundary conditions of the motion are maintained steady. 
Let its velocity and magnetic field intensity be u; and H;. Then 
the time rate of change of dissipation of energy is 


oQ 0 5 
ap Oienberan CAN Rede ooh (3.1) 
D 
We interchange the differentiation and integration and obtain 
0Q fe 
aa Tigui,ge + Tig,eti,g + GB lilist dv. (3.2) 
D 
Replacing primes by ( ),¢ in (2.4), we find 
Ty, Mig = Tigvi,ge. (3.3) 
Equation (3.2) then becomes 
0Q ie 
Paes = 2| (Temes a 14) i,t/0) dv. (3.4) 
| D 
Similar to the verification of (2.6) and (2.15), we find 
0Q 
OE = — 2 | (nuesttas oh iy tH,, 33/0) dv (3.5) 
D 


since u;,4 = 0 and H;,4 = 0 on the boundary. 
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Using the governing equations of magnetohydrodynamics 2) 


pui,e + pujzuig = pPy + Li — fa + 4,43, 


3.6 
Hij.t =F Uj;H 4,5 = ac Te a Ayuy,; =i) 
and (2.21), we have 
eQ l 
ee 2 || pusews,e + cm A;,+Hj,t ) dv. (3.7) 


D 
The integral is non-negative, implying that the dissipation di- 
minishes continuously until @Q/ét = 0. @Q/ét vanishes only when 
uit = O and Hj = 0, ie. when a steady state is reached. Hence 
when the boundary conditions remain steady on &, the motion in D 
tends to its steady state if the motion satisfies (2.21). 

Furthermore, the ultimate state, the steady motion, is unique. 
Suppose an unsteady motion to lead to two different ultimate 
motions M@) and M@). The existence of conditions (2.21) assures 
the linearity of the governing equations. It means that the difference 
of these two motions is also a possible motion with u; = 0 and 
H; = 0 on the boundary. But such a motion necessarily comes to 
rest. The motions M® and M‘) must tend to become and to 
remain identical. If one of these two motions is steady, the other 
must finally come to coincide with it. This assures the uniqueness. 
Hence, the steady motion with equipartition of energy is a stable 
one. In conclusion, we have 

Theorem IIA. In magnetohydrodynamics an unsteady motion of 
an electrically conducting viscous incompressible fluid with boundary 
conditions being maintained steady will always tend to its steady 
state, which is stable and unique, provided conditions (2.21) are 
satisfied. 
_In the absence of a potential Q of the extraneous force we can 
similarly prove the following theorem: 

Theorem IIB. In the case of the non-existence of a potential Q of 
the extraneous force, theorem IIA is true, provided the motion satisfies 
conditions (2.27). 


§ 4. Discussions. In the previous sections we have shown that 
a motion of an electrically conducting incompressible fluid satisfying 
(2.21) or (2.27) has an absolute minimum of energy dissipation. 
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The result is applicable to ordinary hydrodynamics by letting the 
magnetic field be zero. It reduces to the theorems discussed by 
Helmholtz, Korteweg and Rayleigh. The theorems may also 
be used in magnetohydrostatics. In the study of hydrostatic 
equilibrium of magnetic stars, Chandrashekhar and his as- 
sociates in a series of papers ’) 8) have discussed the stability of a 
force-free field, which means that (curl H) x H vanishes. Clearly, 
the force-free field satisfies (2.21) ; it possesses a minimum dissipation 
of energy. 

By the inspection of the governing equations of magnetohydro- 
dynamics, it is seen that when (2.21) or (2.27) hold, the equations 
are linear ones. Superposition is therefore applicable. A discussion 
by Kapur ®) on the possibilities of superposition of solutions has 
led to (2.21) in somewhat less general form. 

When (2.21) or (2.27) hold, a minimum principle is involved. 
This readily suggests that the variational method may be used in 
some magnetohydrodynamic and magnetohydrostatic problems. 
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